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Abstract 

Several  scattering-theoretic  formulations  of  the  theory  of  solids 
are  used  to  discuss  the  energy  levels  associated  with  inperfectlons , 
surface  states,  and  electron  energy  levels  in  highly  disordered  sub- 
stances. The  theoretical  bases  are  thoroughly  developed  and  illustrated. 
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INTRODUCTION 

This  Progress  Report  concerns  the  theory  of  Imperfect  solids,  with 
emphasis  on: 

(I)  atomic  Imperfections  (and  their  optical  properties) , 

(II)  effects  of  proximity  to  a surface  on  these  imperfections  (e.g. 
"recombination  centers") , 

(III)  the  nature  of  surfaces  proper,  and  finally, 

(IV)  transport  in  imperfect  solids. 

It  is  in  the  nature  of  a technical  foray  into  solid  state  theory,  to 
illustrate  current  applications  of  scattering-theoretic  techniques  to 
the  study  of  vacancies,  electron  states  in  highly  disordered  solids, 
the  mathematical  formulation  of  surface  states,  etc.  Other  applica- 
tions: small  polaron,  excitons,  impurity  bands,  etc.  will  be  examined 

in  subsequent  reports. 

Additionally,  a forthcoming  report  will  deal  specifically  with  more 
practical  aspects  of  this  research:  the  study  of  the  photoelectric  de- 

composition of  1^0  ("photolysis")  by  means  of  n-type  TiO^-based  solar 
cells  or  by  similar  other  materials.  The  topics  treated  in  the  present 
report  have  been  inspired  by  the  study  of  these  "more  practical"  aspects. 
Other  immediate  applications:  the  study  of  catalysis  and  corrosion  be- 

gins with  ‘a  model  foreign  atom  or  molecule  or  defect  proximate  to  a sur- 
face; such  models  can  be  anslyzed  by  the  methods  illustrated  in  Part 
III  below,-  as  we  shall  show. 
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I.  ATOMIC  IMPERFECTIONS  IN  AN  INSULATOR 


In  the  Wannier-fimctlon  representation,  a foreign  atom  at  site  Rq 
introduces  a diagonal  perturbation  VQ  as  well  as  a change  in  the  over- 
lap matrix  elements  to  the  nearest -neighbors.  If  the  magnitude  of 
the  former  Is  sufficiently  great,  or  of  the  latter  sufficiently  reduced, 

ic 

a bound-state  is  formed.  If  V >0  the  state  comes  out  of  the  top  of 

o r 

the  valence  band,  the  bound  state  known  as  an  acceptor-level.  If 

V < 0 it  condenses  out  of  the  bottom  of  the  conduction  band  and  creates 
o 

a donor-level.  If  the  point-group  symmetry  of  the  solid  Is  lowered  by 
the  Impurity,  the  conduction  and  valence  bands  may  both  contribute  to 
the  bound  state  level  and  create  a "re combination -center"  or  trap.  The 
symnetry  may  be  lowered  either  spontaneously,  or  by  a Jahn -Teller  dis- 
tortion, or  by  proximity  to  the  surface  (cf.  Chapter  II,  below),  each 
having  different  optical  consequences. 

The  vacancy  may  be  treated  (albeit,  unconventionally)  by  the  same 
methods  as  apply  to  Impurity  atom.  Consider  the  F-center  (a  missing 
Cl  atom  in  NaCl  crystal):  while  conventional  studies  have  created  the 

bound  state  out  of  Na  orbitals  on  atoms  in  the  vicinity  of  the  vacancy, 
we  find  it  simpler  to  postulate  the  continued  existence  of  Wannier  or- 
bitals centered  at  the  position  of  the  vacancy  Rq,  shifted  upward  by 
an  amount  Vq  due  to  the  modification  of  the  Made lung  potential  at  that 

* 

A long-range  potential  such  as  the  Coulomb  potential  has  an  Infinite 
number  of  bound  states;  the  occupancy  being  by  one  or  at  most  two 
electrons,  the  extra  states  are  not  usually  of  crucial  Importance. 

Thus  we  use  the  zero-range  model,  in  which  the  perturbation  is  lo- 
calized at  the  impurity,  in  this  work. 

This  point  has  been  independently  noted  by  others:  cf.  Jaros  and  Brand, 
Phys.  Rev.  B 14,  4494  (1976). 
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site,  and  having  overlap  Integrals  Koj  modified  from  the  values  appro- 
priate to  the  levels  of  an  ordinary  Cl  atom.  The  vacancy  Wannier-orbi- 
tal  will  be  ls-like  rather  than  3p-like  (there  is  no  core  of  occupied 
states  at  the  vacancy,  hence  no  exclusion-principle  requirements). 

The  0 vacancy  in  reduced  Tit^  is  equally  a color  center,  although 
more  complicated  because  it  is  required  to  accomodate  2 electrons.  The 
problem  of  2 Interacting  electrons  orbiting  a comii.on  attractive  well 
has  been  solved  exactly,  and  it  has  been  found  that  if  the  2-body  re- 
pulsive potential  U exceeds  a critical  value  U£  both  electrons  cannot 

be  bound.  But  even  when  U < U , the  first  electron  is  easier  to  ionize 

c 

than  the  second. 

Thus  taking  as  a simple  model  of  an  impurity  or  vacancy  the  one 
first  studied  by  Slater  and  Koster  some  20  years  ago  (VQ  4 0, 
same  as  in  the  absence  of  the  imperfection)  we  have: 


* JV  <,n<r-Rl>  «'<r>  V'V 


i0 

Sn.m  ' \ 


for  i or  j 4 o 
if  i = j ■ o 


because  H',  0^,  0^  are  highly  localized.  Moreover,  if  H' (r)  is  reason- 
ably constant  in  the  atomic  cell  at  R then  g ~ g 6 , is  approxi- 

o n,m  n n ,m 

mately  diagonal  in  the  band  index  (a  constant  potential  doesn't  mix 
atomic  levels) . (We'll  also  consider  the  consequences  of  band  mixing 
which  are  ■very'  interesting).  Let  g^  s coupling  constant  (strength  of 

*D.C.  Mattis  and  E.H.  Lieb,  J.  Math.  Phys . 7,  2045  (1966). 

**G.F.  Koster  and  J.C.  Slater,  Phys.  Rev.  95.  1167  (1954). 
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perturbing  — potential,  previously  denoted  V ) . If  a substituional 
impurity  is  more  attractive  to  electrons  than  a regular  host  atom  that 
it  replaces,  then  g < 0;  if  less  attractive,  g > 0. 

We  can  use  definition  of  Wannier  functions  to  find  matrix  repre- 
sentation in  the  Bloch  states  (nk,mk') : 


^ ^nk,mk'  N ®nm  6 


i(k'-k) .R 


JL  » 

■ — e o e 
N Bn n,m 


We  shall  now  do  the  following: 


i(k'-k) *R 


for  g diagonal.  (1.2) 
n,m 


1.  Obtain  bound  states  of  H + H' 

o 

2.  Obtain  scattering  states  of  HQ  + H' 


3.  Obtain  Scattering  Cross  Section  of  these  ('W>bility") 

4.  Obtain  Optical  matrix  elements  connecting  the  bound  states 
with  any  of  the  conduction  states,  and  obtain  resulting  optical  spec- 


trum, that  is,  do  a complete  analysis  of  "acceptor"  (g  > 0)  or  "donor 
(g  < 0)  levels.  At  the  end,  we'll  consider  gnm  not  diagonal,  for  an 
exactly  solvable  theory  of  recombination  centers  or  traps . 

For  now,  let's  start  with  bound  state  (if  any)  which  takes  the 


*1 


$ (r)  = linear  combination  of  all  Bloch  states 

n,°  in  nth  band 


7n  k ^n.k  ^nk^ 


(1.3) 


We  use  Schrodinger  equation: 


(H  + H')#  - E # 

' o 'no  no  no 


(1.4) 


to  get  the  coefficients  and  eigenvalue  E 


no 


o no  /N  ^ nk 


/N  5 ^nk  Fnk  '•’nk 


+ H*$  ■ g E e 

n0  n3/I  n kk* 


i(k’-k) *R 


Fnk ' ’•'nk ' 


(1.5) 


E $ = tet  E F , E t, 

no  no  /N  nk  no  nk 


As  the  are  orthonormal,  the  coefficients  must  be  equal  on  both  sides 
of  the  equation. 


Fnk(E„o  - Enk>  ■ *n  H S . e 


. i(k*-k) .R 

1 ^ o 


Fnk* 


(1.6) 


So: 


ik.R 


nk  E -€  . n 
no  nk 


where 


. ik'.R 

An=N?,e  Fnk' 

“ ‘ gn  An  N € ,7-E 


nk  no 

= - g,  A S (E  ) 

”n  n n no' 


(1.7) 


which  serves  to  define  Sn(E)  . Either  l 

terminlng  E ) or 
no 


Ar  - 0 (trivial  solution) 


~g_S  (E  ) is  satisfied  (de-^ 
u n v no'  v , 


(1.8) 


Sn (oo)  is  generally  a complex  function  of  its  argunent.  Using 


Im  — * n6  (X)  we  have : 
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— 


and  also 


(1.9) 


Thus  A equation  has  no  solution  for  E in  continuum,  because  r-h-s  is 
n ^ n 

complex.  If  g <0,  there  is  a solution  below  continuum,  1 = -g  R at 
v °n  °n  n 

"A"  provided 


8nl  > rr  (R"in  V R»(W> 

min 


If  ^ > 0 there's  a solution  above  continuum  if 


t>ir 


max 


Analytical  Example: 

To  fix  ideas,  let 

Pn(u>)  " “ Jk  -U)2 , |u>|  < h (1. 10] 


proper  behavior  near  band  edges 
but  omits  van  Hove  singularities 
Bandwidth  * 1 ■ unit  of  energy. 


0)  s % cos  0,  S(<u)  *=  -4e  , 

and  we  verify 

I ■ 4 sin  9 ■ rrp 
n rn 

For  tu  = ± h cosh  X, 

S(u>)  = =F  4e  ^ , with  Im  S = 0. 


The  eigenvalue  equation  becomes: 


1 - (-gn)(T  4 e“^) 
or 

|X|  = ln(±4gn)  for  1 gj  > 


End  of  Example. 


For  the  continuum  solutions  we  take 


*nk  " An,k^nk  + N Lk,k'  *nk'^ 


(1.12) 


where  A^  = normalization  constant. 


Ak-  '1  + \iLk>k.i2ri- 1+  o(i) 

N 


(1.13) 


Equating  coefficients  of  f ^ on  t>°th  sides  of  the  Schrodinger  equation 


yields ; 


E„k  ‘ €„k  + + W~l  Lk,k'> 


(1.14) 


For  many  purposes  = + 0(1/N)  can  be  replaced  by  — i.e.  the 

continuum  of  exact  scattering  solutions  "interlaces"  the  continuum  of 
unperturbed  Bloch  states.  Next,  equate  coefficients  of  : 

ik-R  ik".R  -ik'-R 

lk,k'(Enk-«„k'>  * *„<e  0 <r-15> 

A series  of  obvious  steps  follows  (solve  for  in  terms  of 

-1  ik".R0 

D ■ N k"  6 then  solve  for  D self -consistently  using  the 

value  of  obtained  previously)  easily  leading  to: 


i(k-k')  -R 


lk-k'  ‘ €„k  - <W  1 + 8ns„<enk) 


(1.16) 


Then,  (14)  becomes: 


f”k  + N 1 + e„sn(€nk) 


(1.17) 


We  now  verify  that  the  bound  state  is  orthogonal  to  the  scattering 
states,  and  that  the  latter  are  orthogonal  to  one  another. 


(a)  Orthogonality  of  $nQ  to  $nk: 


The  Integral  of  (1.3)  X (1.13)  yields:  (recalling  are 


orthonormal) : 


lk-R 


lk'.R  i(k-k') »R 
g o o 

n „ e e 


VEo  "N  k-  Sk'-Eo  6k'  -«k  1+E»V6nk> 


lk-R 


VE„  11  ‘ l + 8„S„«»k> 


o 

ik*R 


° ,1+  *nW 


€. -E 
k o 


^1  + 


*nSn^nk> 


This  vanishes  by  virtue  of  the  eigenvalue  equations  (1.7),  (1.8)  which 
determine  E 


no 


(b)  Orthogonality  of  one  scattering  state  with  another: 

i(k-k') *R 


Lk,k'  + L k’k  + N k"  Lk'k,,Lkk"  = 6 


{ — l r 3a 


*n 


«nSn 


nk' 


*h8n  (€nk-> 


] 


*n 


1 + *nSn<enk>  1 + ^V^nk 


T • N z — h — >*(P — IT-)!  - 0 

N k"  €nk"^nk'  ^nk'^nk 


assuming  appropriate  infinitesimal  inaginary  parts  in  each  denominator 


(according  to  the  prescription:  + 10+,  6^,  - i0+  for  k,  and 


similarly  for  the  others,  satisfying  the  requirements  of  causality 
thus) . 

Scattering  by  Impurity 


2 

The  differential  scattering  cross-section  a.  (0)  = | f.  | , where 


-1 


f.  (0)  - lim  rN  E L.  . , ♦ , (r) 

•»  j CO  • •»  | •»  K 


9 . 


(1.18) 


in  which  cos  9 *»  k*k7kk'.  We  neglect  the  periodic  component  of  the 
Bloch  functions,  and  perform  the  sum  over  the  phase  factors  only.  After 


integrating  over  angles  and  approximating  by  k we  have  left: 
f S o.i  1 1 • 2 1 sin  k r 

fk<8>  ~ i+  |s(ek)  rdk  k ^27? 


(1.19) 


which  yields: 


fk(6)  ~ 1 + gS(fk)  6 


(1.20) 


As  this  is  independent  of  0 the  total  scattering  cross-section  is: 


a.  = const,  x 5- 

(i+gR)z+  (g"Pr 


(1.21) 


The  scattering  time  is  related  to  the  scattering  cross-section  by 

and  the  mobility  (the  conductivity  per  carrier  at  a given 
energy)  is  e T^/m  , where  v^  = and  1/m  = (suitably  aver- 

•ff 

aged).  For  a parabolic  band  approximation,  v^  « p(f^)  and  m = const. 
This  is  sufficient  for  our  purposes.  For  N^.  impurities  (concentration 
Cj.  ■ Nj/N)  we  have : 


u 


2n  g p(ek) 

1 (1+  gR(€n))2  + (gTTP(€k))2 


(1.22) 


after  juggling  the  multiplicative  constants  to  obtain  agreement  with 
the  first  Bom  approximation  (Fermi's  Golden  Rule)  at  g = 0.  An  iden- 
tical result  is  obtained  directly  from  the  relation: 


2 r ■ »I 


(1.23) 
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with  Enk  given  in  (1.17).  The  general  dependence  of  scattering  time  on 
energy  is  Indicated  in  the  figure  below.  Several  remarks  are  also  in 
order: 

(1)  the  Born  approximation  is  very  bad  when  there  exists  a bound 
state,  even  at  the  band  edge,  where  p -•  0.  The  denominator,  instead 
of  being  1 is  (1  + gR)^  = g^(R(^k) -R(Eq))^  which  can  be  either  much 
larger  or  much  smaller  than  1;  for  a bound  state, Eq,  very  close  to  the 
band  edge  (as  for  a typical  donor  or  acceptor  level),  this  expression 
will  tend  to  be  very  small  at  the  band  edge,  resulting  in  true  scatter- 
ing being  much  enhanced  over  the  approximate  Bom  value. 

(2)  At  the  resonance  energy  (^kres  " ^min  « ^min  “ E0^  t*ie  term 
2 

(1  + gR)  in  the  denominator  vanishes  entirely,  and  the  scattering  time 

2 

becomes  extremely  short;  the  coupling  constant  g in  the  denominator 
2 

and  the  factor  g in  the  numerator  cancel.  We  have  there  (1/t) 

max 

(3)  When  the  potential  is  repulsive  so  that  there  is  no  bound  state 

below  the  Conduction  band  and  no  resonance  at  low  energies,  the  effec- 

2 2 

tive  scattering  strength  is  g /(I  + gR)  at  the  lower  band  edge,  which 

2 

may  be  considerably  reduced  from  the  Born  value  g . Thus  a repulsive 
potential  which  may  be  strong  enough  to  create  an  acceptor  level  near 
the  valence  band  may  Induce  very  little  scattering  in  the  conduction 
band . 

These  remarks  summarize  some  important  aspects  of  scattering 
theory,  which  can  be  explicitly  verified  for  our  soluble  model. 


Some  Optical  Properties: 

The  energetics  for  an  insulator  + 1 impurity  is  the  following: 


assuming  an  attractive  potential  (g  < 0).  The  level  can  be  occupied 
or  empty:  thus  the  impurity  introduces  two  new  processes  for  the  ab- 

sorption of  light,  labeled  A and  B: 


H - Eo 

h>  . > E + E 
A ~ g o 


3 


For  light  polarized  in  the  z-direction  the  optical  matrix  element  is: 


H , - - rr  (A  eiu,t  + A*  e”iu)t)  p , A - — 6 (1.2- 

e 1 -opt  me  u)  u>  z u>  u>  u> 


This  weak  coupling  Hamiltonian  is  currectly  treated  in  lowest  Born  ap 
proximation.  For  long  wavelength  light  (not  X-rays)  we  may  assume  A 
to  be  constant  in  space,  and  obtain  for  the  probability  per  unit  time 
of  process  A or  B occuring: 


u> 


, |£  |2c2 

" * *“  <=»  (-“-r-)  ,,Ej<fin.i|p,|i„itt.i)|2  «<eflnil  ± 


0) 


final 


ft  <u-f 


final5 

(1.25) 


< 

i 
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with  (+)  for  emission  of  photon  and  (-)  for  absorption.  For  the  sake 
of  definiteness  we  treat  the  absorption  spectrum  only,  the  reader  can 
make  the  necessary  modifications  to  obtain  the  emission  spectrum.  Be- 
cause the  spatial  dependence  of  the  electromagnetic  field  is  Ignored, 

we  require  only  the  matrix  elements  of  p , which  are: 

z 


<l>z>„'k\„k  E JV  +„'k'<r»>z  Wr> 


’ Vk'Mn'n(k> 

(Note:  conservation  of  k) 


(1.26) 


For  n’=n,  we  have  the  simple  result: 

M ■=  mv  (n,k)  (1.27) 

n ,n  z 

where  v is  the  z -component  of  the  velocity: 

-K€„,k  v-2* 

For  n'  4 n there  exists  an  f-sum  rule  relating  the  matrix  elements  of 

p to  one  compoent  (zz -component)  of  the  effective  mass  tensor: 
z 


E 

nVn 


(1.29) 


where 


£ ' ^Pz^n'k,nk' 

m € i,  * € . 
n k nk 


(1.30) 


A.H.  Wilson,  "The  Theory  of  Metals",  Cambridge  University  Press,  p.  46, 
1954. 
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The  other  components  and  a proof  are  given  in  Wilson,  op.  cit.  In  an 
approximate  model  where  only  2 bands  are  considered,  this  sum  rule  re- 
lates the  curvature  of  the  bands  to  the  optical  matrix  element;  so  the 
parameters  in  the  model  are  closely  related  and  are  not  to  be  chosen 
arbitrarily. 

Now,  either  Mn,  ^(k)  Is  zero  at  k = 0 (forbidden  transition,  e.g. 
from  an  s-like  band  co  a d-like  band,  both  of  which  have  the  same  par- 
ity), but  it  won't  be  zero  at  k 4 0:  thus  M , (k)  = constant  x k is  a 

n'n  z 

good  estimate  for  its  value  at  k = 0.  Or,  M , (0)  is  non-zero  at 

n ,n 

k = 6,  as  for  transitions  between  an  s-like  and  a p-like  band.  In  the 

case  of  such  allowed  transitions,  we  shall  approximate  Mnin00  “ 

M , (0)  by  its  vale  at  k = 0,  for  simplicity,  and  treat  the  matrix 

n ) n 

element  as  a constant. 

In  the  study  of  optical  properties,  it  is  frequently  helpful  to 
use  the  identity 


• iJV  C'k  8 *nkl 


IJV  ♦i'k  8 "I  “I 


(1.31) 


where  the  last  result  uses  energy  conservation,  and  is  therefore  only 
correct  for  use  in  equation  (1.25) . 

We  can  now  study  process  A.  Label  the  initial  states  (v)  for  "val- 
ence band",  the  final  states  (c)  for  "conduction  band".  They  are,  re- 


spectively: 
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‘vk-  *vk+  N l<V)k,k’  *vk' 


where  Lv 


i(k-k') -R 

» . e 1 8v 

k,k'  €vk  " €vk'  1-gvSv((W 


(Note:  see  (1.16)) 


(1.32) 


where 


'c,o  ■ TsT  2 Fck  *ck 


-ik*R  g A 

Fck  = e ° » (by  (I,7) 

o ck 


,8sC(Eo\-% 

*c4c  ■ <— §E > 

o 


Thus,  the  matrix  element  is: 


(1.33) 


-Ik*R  , , 

<*vklPj*co>  = 6 ° X 


c o o 


M (k)  . . g.  M (k*) 

i_vc i r , _ i k* f*  cv  , 

Vck  Bk'«vk^vk'  i-«,SiKvk>  Eo^vk' 


(1.34) 


If  the  transition  is  forbidden,  M (k’)  is  odd  in  k'  and  the  second 

cv  z 

term  in  curly  brackets  (the  sum)  will  vanish  by  symmetry.  Even  if  the 
transitions  are  allowed,  for  g < 0 (as  we  are  assuming)  there  is  no 


resonance  near  the  top  of  the  valence  band,  so  the  second  term  contri- 
butes little  to  the  absorption  near  threshold.  We  drop  it  for  simplic- 
ity, obtaining  for  the  transition  rate  per  impurity  atom: 


«-r<5?2 

i M (k)  - 
1 „ , cv'  '2 


2 x OS  (€  )OE  )_1  x 
' c o o' 


»i(M  s(Eo^.k> 


(1.35) 
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r 1 
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X 

which  is  readily  seen  to  behave  as  (a>-u>T)  near  threshold  when  the 


.3/2 


transitions  are  allowed,  and  (cu-u^,)  when  forbidden,  rising  to  a peak 

-4 


before  falling  off  as  <u 


Before  this  happens,  absorption  at  the 


fundamental  absorption  edge  fiw  - E has  already  obscured  the  line. 

S S 


From  (1.25)  and  (1.26)  the  fundamental  absorption  is  given  by 


w = r O Kr  = Kv<k>l  *(erk  - - evk> 


U) 


cv 


ck 


"vk 


(1.36) 


It  starts  at  E , also  rises  as  (u)-u>  )^  if  M is  constant,  or  as 
g g cv 

3/2 

(<B-cu  ) if  the  transition  is  forbidden,  but  the  overall  absorption 
8 


rate  is  a factor  N/Nj.  larger  than  the  total  impurity  absorption  (every 
site  in  the  crystal  contributes  to  the  fundamental  edge,  whereas  only 


the  Nj.  impurities  contribute  to  the  impurity  process;  thus  the  impurity 


is  visible  only  in  the  range  of  frequencies  where  the  pure  crystal  is 
otherwise  transparent) . 

We  next  turn  to  process  B;  if  the  bound  level  is  close  to  the  con- 
duction band,  the  optical  line  associated  with  this  process  will  be 
reasonably  narrow  and  quite  low  in  frequency  as  compared  to  the  funda- 
mental absorption,  thus  this  process  will  be  more  distinguishable  from 
the  background. 


The  absorption  matrix  element  is:  MRn  = mvz  (equation  (1.4)),  so 


the  absorption  rate  per  Impurity  becomes: 


m2v  2(k) 

r <fc>2 1*/  >«ck^-E0> 

(E  -€  , ) 
o ck 


The  sunned  term  in  the  matrix  element,  E L.  .|V  (k')  (E  -€  . ,) 

1.1  «Cj  K Z O ClC 


-l 
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3/2 

also  satisfying  a (<j)-(0o)  law  at  threshold.  The  various  results  are 
plotted  in  the  figure  on  page  19. 

We  note  that  while  all  the  formulas  in  this  section  have  been  ob- 
tained without  any  significant  approximation,  nevertheless  the  model 


is  not  accurate.  At  the  least  we  should  include  the  modification  in 
overlap  integrals  discussed  on  page  3,  as  well  as  2-electron  repulsion 
to  have  a reasonable  model  with  which  to  compare  optieal  properties  of 

0 vacancies  in  Ti02  for  example.  As  this  requires  no  fundamental  de- 
parture from  the  methods  of  computation  given  here,  we  shall  no  expand 
on.  this  but  shall  merely  give  the  results  in  a subsequent  report. 


i 


II.  EFFECTS  OF  SURFACE  PROXIMITY 


I 


Proximity  to  a surface  will  be  seen  to  affect  the  binding  energy, 
to  cause  structure  in  the  optical  absorption  spectrum,  and  to  lower  the 
point-group  symmetry  of  the  impurity  site  to  where  it  can  act  as  a re- 
combination center. 

For  the  purpose  of  this  chapter,  we  shall  assume  that  the  boundary 
condition  typifying  a surface  is  \)i  = 0 at  z = 0.  The  solid  will  con- 
sist of  atomic  planes  at  z = ns0>  n=l,2,3,...  . An  impurity  at  a depth 
d = paQ  from  the  surface  can  be  rendered  soluble  by  the  artifice  of  an 
"image"  impurity  at  -pa^.  The  eigenstates,  whether  bound-state  or 
scattering  states,  are  then  computed  for  the  infinite  solid  (-»  < n < +“) 
but  only  functions  odd  under  reflection  about  z = 0 are  retained,  thus 
automatically  satisfying  the  boundary  condition. 

We  examine  the  bound  states  of  the  impurity  and  its  image,  follow- 
ing the  procedures  of  chapter  I.  Schrodinger 's  equation  assumes  the 
form: 


H$ 


;k*k(r> 


t£y(k'-k)-d  + 


-i)k'-k).d 


’ -4<d>  75  E Vk  <n-l> 

in  which  A(d)  is  the  distance -dependent  binding  energy  (i.e.  A(°>)  = 
|Ej)  . Equating  coefficients  of  yields: 


Fk<A  + €k)  » jj-E,  FR|  2 cos(k'-k)  -d  = g(e"ik'd  ?(d)  + c.c.)  (II.2) 


2t 


in  which  we  define 


1 „ „ ik'.d 

s E Fk-  e 


(It. 3) 


The  self-consistent  solution  of  (II. 2)  and  (II. 3)  yields: 


|(±d)  = g[(§(±d)T(0)  + ?(=Fd)T(±2d))] 


(II .4) 


where 


. ik’ «x 

T(X)  H N ^ A(d)+€k, 


(II. 5) 


The  pair  of  equations  (II. 4)  can  have  a solution  only  if  the  secular 
determinant  vanishes: 


l-gT(O) 

-8T(-2d) 


- gT(2d) 
l-gT(O) 


(II. 6) 


which  is  equiva lent  to 


1 - g(T (0)  ± T(2d)) 


(II. 7) 


Only  the  solution  belonging  to  (-)  is  physically  acceptable,  as  it 
alone  satisfies  the  boundary  condition;  it  corresponds  to  ^(d)-  -£(-d) 
and  thus  = 0 at  d * 0 according  to  equation  (II. 2).  The  (+)  solu- 
tion is  even  at  about  the  plane  z » 0,  and  must  be  discarded.  Thus  the 
eigenvalue  equation  reduces  to: 


1 1 _ l-e21k>d 

g " N l A (d)+f  k 


(II. 8) 


We  shall  now  establish  that  it  is  entirely  possible  for  g to  satisfy  the 
condition  for  the  existence  of  a bound  state  in  the  bulk,  yet  fail  to 


satisfy  (II. 8)  near  the  surface. 


22 


For  this  purpose,  it  is  sufficient  to  calculate  gc(d)>  the  minimum 

g which  will  yield  a bound  state  according  to  equation  (II. 8)  above. 

Deep  in  the  bulk,  at  d » ®,  upon  setting  A - 0 we  have: 

o 


_1_  . I S 1.  . s(0)  . R(0) 


(II. 9) 


with  S(id)  previously  defined  in  chapter  I.  The  change  upon  nearing  the 
surface,  formed  by  setting  A(d)  **  0 in  (II. 8),  is: 


i i i 2ik-d  o 

__i I . I y £ = 2m 

g (»)  g (d)  N L € 4^2df 


(II. 10) 


in  which  m = effective  mass  of  electrons,  and  the  integrand  is  recog- 
nized as  the  Fourier  transform  of  the  Coulomb  potential.  Solving  this 


equation  for  gc(d),  we  have: 


gc(“) 


8c  ^ l-(mg  (»)/4rrjd 


(II.  U) 


The  results  are  shown  in  the  figures.  The  binding  energy  decreases  as 
the  surface  is  neared,  and  may  even  vanish  for  Impurities  located  within 
the  first  several  atomic  layers. 

The  effects  of  proximity  to  a surface  on  the  optical  absorption 
spectrum  are  even  more  dramatic.  We  obtain  them  by  first  solving  for 


, ,v  2i  sin  k«d 

Fk  - Tw  V fk 


(11.12) 


Normalization  requires: 


1-h  |f  I2  - 2g2|S(d)|2  i E -1-"-03  ?k'd-  2 

N k N (€k+A(d))Z 


(11.13) 
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•er  i*J(* 


Plot  of  gc(n)/gc(®)  as  schematically  obtained  from  equation  (II. 11). 
For  a typical  attractive  potential  g a bound  state  may  be  produced  in 
the  bulk  but  not  in  the  surface;  for  the  illustrated  value,  the  criter- 
ion for  a bound  state  is  met  only  for  the  surface  layers  n > 3,  no  bound 
states  exist  if  the  impurity  is  within  either  of  the  first  two  atomic 


2 


By  use  of  (II. 8)  this  can  be  written  as: 


Uc«i)l2  - %82  Sff- 


(11.14) 


and  can  be  used  to  calculate  the  absorption  matrix  element: 


|(t  It  )|2  = Lh2 

1 wk'V’  N dd  (A  + €.) 


(11.15) 


thus  yielding  for  the  optical  absorption: 


,,  sin  2kd* 
“ ’ "»<1  - 2kd  ^ 


(11.16) 


in  which  w^  is  the  transition  probability/unit  time  given  in  chapter  I, 
with  A (d)  replacing  A (“) . The  extra  structure  comes  from  the  trigono- 
metric term,  which  in  the  effective -mass  approximation  has  the  value 


k = (2m)' 


(a)-Eg  + 


A (d)  ) 


(11.17) 


Band  Mixing: 

A second  important  effect  of  the  proximity  of  a surface  relates  to 
the  lowering  of  symmetry.  Where  a symmetric  perturbing  potential  would 
not  mix  valence  and  conduction  bands  (particularly  if  these  have  in- 
trinsically opposite  parity)  the  impurity  '‘molecule"  (imperfection  + 
mirror  imperfection)  near  the  surface  generally  has  a lower  symmetry. 
Valence-  and  conduct ion -band  states  are  mixed  into  the  bound  level, 

m 

which  now  becomes  equally  accessible  to  electrons  and  to  holes  — i.e. 
a recombination  center. 

A model  for  this  starts  with  the  following  perturbing  Hamiltonian: 


HVk-,„k  - I + V,n<1-5n,„'>l8ln  k'R.  sl"  V (II’18> 
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Restrict  n to  2 bands  "c"  and  "v"  and  obtain  the  bound  state  if  any. 
The  bound  state  wavefunction  must  take  the  form: 


§ = N ^ E F .ill  . (r-R  ) 

o nk  n,k  n,kv  o' 


(11.19) 


in  which  \|i  , (r)  is  the  linear  combination  of  two  Bloch  functions, 

n,K 

(k  ,k  ,k  ) minus  that  at  (k  ,k  ,-k  ),  which  vanishes  on  the  z “0  plane, 
x y z x y z 

Schrodinger 's  equation  now  reads: 

HI  = T:  S F , € , ♦ , + Tv  Z F , € . t, 
o /N  k vk  vk  vk  /N  ck  ck  ck 

2 8 2g 

+ E sin  k.R  sin  k'*R  F .ill  , , + --y /»-  E sin  k.R  sin  k'*R  F ,i|i  . 

N3/2  kk'  0 o vk  vk'  n3/2  kk,  o o vk’ck 


2gc  2gvc 
+ rz  E sin  k.R  sin  k'‘R  F ,f  £ sin  k*R  sin  k'  *R  F 

n3/2  kk.  o o ck  ck'  n3/2  kk,  o o 


(11.20) 


and  we  equate  coefficients  of  iji^  and  iji^^  in  turn.  For  the  first: 


(Eo^ck)Fvk  " 2 Sin  k'R~[8„  h S.F.^.sin  k’*R.  + s,.„  k k'.RJ 


o L °v  N u vk 


o vc  N 1,1  ck 


s 2 sin  k*R  [g  D + g D ] 

O ^ V V VC  c 


(11.21) 


which  serves  to  define  D . For  the  second: 

n 


<Eo-eck)Fck  ■ 2 sl"  k-Ro  [8cvDv  + SA1  (n-22> 


Evaluating  D using  the  definitions  above: 
n 


D •=  jr  E (E  -f  . ,)"1i(2  sin2k'  *R  ) (g  D +g  D ) j = -S  (E  ) [g  D +g  D ] 
V N kt  o ck  1 O V V VC  C 1 V O V V vc  CJ 


_ 2 sin  k.R 

VE>  - »'  l -€-(»  (II‘24> 

by  analogy  with  (1.7).  Similarly  for  Dc: 

D = -S  (E  )[g  D + g D 1 (11.25) 

c c o c c cv  v ' 

These  homogeneous  equations  (11.23)  and  (11.25)  have  a solution  iff  a 
secular  determinant  vanishes.  A special  case  of  this  is  trivially  sol- 
uble: let  all  g's  be  equal  (i.e.  strong  interband  mixing),  then  the 

condition  for  the  bound  state  (recombinatior  center)  to  exist  within 
the  energy  gap  reduces  to  the  simple  algebraic  form: 

g"1  = -[Sv(Eq)  + Sc(Eo)]  (11.26) 


We  plot  the  r.h.s.  of  this  equation  as  well  as  the  d.o.s.  functions  of 
the  two  bands  in  the  Figure: 


The  bound  state  evidently  has  a dipole  moment  if  the  two  bands 
have  opposite  parity.  The  magnitude  of  the  dipole  moment  can  be  ob- 
tained through  (1.31)  as  a function  of  Mvc(k) , We  omit  this  calcula- 
tion, but  note:  a first-order  Stark  effect  is  implied.  That  is,  the 
bound  state  will  respond  linearly  to  an  applied  electric  field  fi  while 
the  bulk  state,  in  which  the  bands  are  not  mixed,  is  affected  only  to 
0<?2)  in  such  a field.  This  should  help  to  identify  recombination  cen- 
ters optically,  by  the  techniques  of  electro-reflectance. 


III.  CONCERNING  ELECTRON  STATES  NEAR  SURFACES 


[ 


The  study  of  surfaces  has  become  a very  important  and  active 
branch  of  solid  state  physics.  Surface  physics  is  particularly  impor- 
tant in  the  understanding  of  solid  state  devices,  including  solar  cells, 
also  for  the  study  of  heterogeneous  catalysis,  including  photolysis,  all 
for  a variety  of  reasons.  The  methods  by  which  surfaces  have  been 
studied  are  typically:  numerical  computation  of  a desired  quantity 

(e.g.  band  structure)  for  a thin  slab  of  N atomic  planes,  varying  N and 
expanding  the  property  as  XN°  + YN^  . Clearly,  X is  the  surface-related 
property.  Just  as  clearly,  this  procedure  is  tedious  and  does  not  lend 
itself  to  rapid  insights.  We  illustrate  in  this  chapter  a scattering- 
theoretic  method  designed  to  focus  specifically  on  the  surface -re  la  ted 
properties  in  a solid  where  N is  immediately  taken  to  be  large  or  in- 
finite. By  way  of  illustration  we  shall  determine,  under  what  condi- 
tions will  electronic  surface  state  energy  bands  be  created. 

Our  model  problem  starts  as  follows:  we  find  the  linear  com- 

bination of  Bloch  functions  to  satisfy  the  boundary  conditions  \|i  = 0 at 
z < 0,  and  evidently,  the  Schrodinger  equation  within  the  solid.  This 
analysis  is  extremely  illuminating  and  has,  apparently,  not  been  done 
before.  We  find  it  Important  to  analyze  the  amplitude  of  the  eigen- 
functions near  the  surface,  especially  near  the  energy  minima  and  maxima 
of  the  band  structure.  We  then  apply  a perturbation:  in  the  simplest 

case,  only  on  the  first  atomic  plane  (l.e.  the  surface  plane).  If  we 
ask  the  following  question:  under  what  conditions  will  an  Infinitesimal 
perturbation  (g  « either  E or  width  of  energy  bands)  produce  surface 
states?,  we  obtain  a remarkably  simple  answer:  only  when,  in  the  2D 
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Brlllouin  zone  of  the  surface,  there  exist  points  or  lines  along  which 


£ 

l:,i 


a certain  component  of  the  inverse-effective -mass  tensor,  or  , vanishes. 

zz 

When  the  surface  perturbation  is  large,  there  always  exists  sur- 
face states,  and  our  method  permits  a relatively  simple  evaluation 
thereof.  In  any  case,  having  obtained  the  surface  bound  states  and 
the  scattering  bulk  states  of  the  terminated  solid,  we  can  proceed  to 
study  the  effects  of  an  added  perturbation  at  the  surface.  (This  per- 
turbation could  be  done  to  a molecule  undergoing  heterogeneous  cata- 
lysis or  photolysis.)  We  shall  see  that  the  effects  of  such  a pertur- 
bation are  qualitatively  different  than  if,  as  in  the  usual  analysis, 
the  bulk  wave  functions  and  band  structures  were  used. 

In  the  periodic  solid  the  complete  orthonormal  set  of  Wannier 
functions  is  related  to  the  complete  orthonormal  set  of  Bloch  functions 
(eigenfunctions  of  the  unperturbed,  periodic  Hq)  as  follows: 


ik-R , 


_ jl  a.K.'K,  , -ik*R, 

Wr>  = N S e J 0n(r-Rj)  and  0n(r-Rj>  = N*.Z  e \fc(r) 

(Bloch)  J (Wannier 

(III.l) 

Thus  they  are  Fourier  transforms  of  one  another  over  the  N discrete 
values  of  R^,  or  k (in  the  space  or  reciprocal  lattice,  respectively). 

The  introduction  of  one  or  more  surfaces  and/or  perturbations  has  the 
effect  of  replacing  exp  ± ik-R^  by  the  components  Uk(Rj)  of  a unitar; 
transformation.  As  we  intend  to  retain  the  as  the  basis,  it  is  im- 
portant to  reformulate  the  operator  Ho»  initially  given  as 

Hq  = -(^/2m)V^  + Vper£0<jic(r)  (with  b = 1 henceforth). 

One  readily  obtains; 


i 
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fd^r  0n*(r-Rj)HQf(r)  = Pn(-id/dR.)F(Rj) 


(III. 2) 


where  F(R^)  *s  the  coefficient  in  the  expansion  of  f,  i.e.: 


f (r)  = E F(Rj)0n(r-Rj) 


(III. 3) 


Equation  (III. 2)  is  valid  for  arbitrary  f(r)  (not  necessarily  an  eigen- 
function of  Hq  or  of  any  other  Hamiltonian) . 

To  introduce  diagonal  (in  the  Wannier-representation)  potentials 

as  perturbations,  we  merely  add:  E V 6^  £ to  H and  finally  ob- 

m m KjjjjKj  o 

tain: 


H Dk(Rj)  - EkUk(RJ) 


(III.*) 


Si 

t! 


as  our  new,  pseudo-Schrodinger  equation,  subject  to  the  boundary  condi- 
tions U. (R.)  = 0 for  Z.  < 0 (for  the  study  of  surface  phenomena;  note 
K J J 

that  periodic  b.c.  could  be  used  if  we  were  not  concerned  with  surfaces; 
indeed,  we  shall  use  p.b.c.  for  the  coordinates  and  Yj) » in  which 
is  the  new  energy  eigenvalue  (which  may  or  may  not  interlace  the 
old)  k is  a quantum  number  which  adiabatically  reduces  to  the  crystal 
omentum  when  H -•  Hq  and  p.b.c.,  are  restored, Uk(R^)  is  the  "wave function" 
defined  on  the  N points  R^  only,  and  finally: 


H - C (-id/dR.)  + E V 6^  - 
nv  y m m R ,R. 

J j 


(HI. 5) 


Is  the  Hamiltonian,  a matrix  defined  'only  on  the  N lattice  points.  It 
should  be  emphasized  that  despite  the  explicit  differential  operator, 
this  H incorporates  only  discrete  translations  by  integer  multiples 
of  primitive  translation  vectors. 
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Surface -Theoretic  Hamiltonian: 

For  surfaces  which  retain  the  bulk  translation  vectors  in  the  X-Y 
plane  the  above  simplify  somewhat: 

H U^Zj)  = EjU^Zj)  (III. 6) 

in  which 

" " WV  •13/8Zj>  + S V„  SZ  ,Z. 

J m’  j 

and  Zm  = raa,  m=l,2,3,...  Note  that  we  have  a different  "linear  chain" 
1/3 

having  N = N distinct  eigenvectors  and  eigenfunctions,  for  each  of 
z 

2/3 

the  N values  of  the  2D  vector  (k  ,k  ) . This  vector  we  denote  kn 

* y II 

hence  forth. 

The  linear  chain  problem  thus  defined  is  far  from  trivial,  as  the 
following  example  will  indicate. 

Example : 

Suppose  the  potential  perturbations  V to  be  absent,  and  consider 

m 

the  solutions  of  (III. 6 and  III. 7)  when  the  band  structure  takes  the 
form: 

€ (kii,k  ) = -2K  cos  k a - 2L  cos  2k  a (III. 8) 

II  II  z z z 

in  which  K and  L are  presumably  functions  of  kjj . A plot  of  ^n(^z)  Is 
given  belo»(kii  is  fixed;  only  positive  k^  are  shown  as  C-  is  symmetric). 


r 

i 

' 


j 

' 


32 


We  wish  to  solve  the  eigenvalue  equation: 


H U(z)  ■ {-2K  cosh  ad/3z  - 2L  cosh  2a3/dz}u(z)  = E U(z)  (III. 9) 

on  the  lattice  points  z ■ ma , m=l,2,...  subject  to  the  b.c.  U = 0 for 
m ■ 0,-1,...  and  U = finite  as  m -*  +®.  We  can  construct  such  a solu- 
tion using  the  set  of  k's  which  satisfy  € (k  ) = E;  label  them  k , 

z q 

q = 1,2,3 ,4.  We  then  have: 


U = E A 


(III. 10) 


and  need  only  determine  the  A^  which  satisfy  the  b.c.  The  solutions 
have  different  character  in  cases  (a) -(c)  pictured  above.  We  start 
with: 

Case  a; 

Fc--.  E in  the  range  -2(K+  L)  < E < + 2(K  - L)  there  are  2 real 
solutions  ± ko  and  2 complex  roots  ± (k^  + ik^) . One  of  the  latter 
grows  exponentially  as  m -•  +®  and  must  be  discarded.  We  find,  after 


some  algebra: 


. -lj-K  , L r .K  . 2 2L'E,%. 

v ■ cos  te + nr r w + i 


»<ki + lk2>  • rra-+ 1 + r<fe>  + 


K .2  . 2L-E  nk 


+ td£i  +«!r>2  + i 


Thus : 


ik  ma  -ik  ma  ik.ma-k„ma 
U(m)  **  Ae  ° + Be  ° + Ce  1 
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satisfies  the  difference  equation  (III. 9)  everywhere  except  at  m = 2 
and  1.  However,  imposing  U'(0)  = U(-l)  = 0 enables  us  to  satisfy  the 
equation  there  also.  The  vanishing  condition  at  m = 0 is; 


A + B + C = 0 


while  that  at  m = -1  requires: 


-ik  a +ik  a -i(k..  + ik  ) 

Ae  ° + Be  ° + Ce  1 2 ‘ = 0 


Normalization  of  the  "wavefunction"  U for  a chain  terminating  at  N 


yields 


|A|2  + 1 B| 2 = 1/N  (in  the  limit  N » 1/k  a) 
z z z 


These  uniquely  determine  a single  solution  for  each  value  of  E.  The 
reader  may  work  out  the  elementary  algebra  for  the  coefficients  A,  B, 
C and  also  verify  that  when  L = 0 this  solution  properly  reduces  to 

k -1 

V = (2/Nz)^  sin  kQam,  where  kQa  = cos  (-E/2K) . 


We  next  turn  to: 


Case  b: 


This  case  is  "non-standard",  in  the  sense  that  the  energy  minimum 
occurs  not  at  k = 0.  In  the  range  € < E < £ the  analysis  is  exactly 

———————  o — — TT 

that  of  case  (a),  and  a unique  solution  with  3 components  is  again 

found.  However,  for  E in  the  range  f < E < 6 , which  encompasses  the 

m ~ — o 

neighborhood  of  the  energy  minimum,  there  are  found  A real  wavevectors 
± kQ  and  ± k^  at  each  energy.  Imposition  of  the  two  b.c.  at  m ° 0 and 
-1  reduces  this  number  to  2 (except  precisely  at  E = where  only  a 
single  solution  survives  — howevew,  this  is  a minor  technical  point). 
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Let  €(kQ)  = €(k^)  = E (cf.  Figure  b) . As  the  k's  are  all  real,  we  use 


sin  and  cos  functions  instead  of  exponentials,  and  guess: 


U^(m)  = ~ | sin  koma+a  sin  k^ma  + ( 


s in  k a + o'  s in  k,  a 
o 


cos  k a - cos  .v, 
o 1 


■)  (cos  k ma  - cos  k,ma)  1 

k,  a'  o 1 ’ 


as  the  first  solution,  and  (m)  = same,  with  k^  and  k^  permuted.  Note 
that  sin(kma)  vanishes  at  m = 0 and  so  does  cos(kQma)  - cos(k^ma).  The 
coefficients  have,  furthermore,  been  chosen  to  satisfy  the  second  b.c. 


U = 0 at  m = -1.  It  follows  that  also  satisfies  both  b.c.'s  It 


remains  to  obtain  the  normalization  constants 


N 


sinka+orsink.a 

C - [<j»  |1  + a2  + 2<^-g-  — Jr1.-  )2t}% 

o 1 


with  obtained  by  the  permutation  of  kQ  and  k^.  Finally,  the  para- 


meter a is  chosen  to  make  orthogonal  to  Uj*  We  obtain: 


Of  = - [2  sin  kQa  sin  k^a]  * x 

2 2 2 2 2 % 
[sin  k a4-sin  k,a4-(cos  k a -cos  k,a)  + j[sin  k a -sin  k,  al^+fcos  k s-cos  k,  a]  t x 
O 1 O I 1 o 1 o 1 1 


2 L 

j[sin  kQafsin  k^a]  + [cos  k^a-cos  k^a]  }"2] 


Finally,  we  turn  to: 

Case  c; 

Treated  as  (a)  in  the  range  £ < E < € . From  £ to  £ we  follow 

O — — TT  TT  m 

(b). 

It  is  now  important  to  examine  the  solutions  near  the  energy  minima 
and  maxima.  We  shall  illustrate  only  the  minima:  the  maxima  are  treated 
analogously,  so  their  analysis  would  be  redundant. 
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. 

We  consider  the  plane  index  m to  be  finite,  and  allow  the  energy 
to  approach  the  minimum  in  case  (a)  to  obtain: 

(a)  U -*  i(2/Nz)^  sin  kma  (as  k -•  0) 

It  is  important  to  note  that  |u|  vanishes  proportional  to  k as  k -•  0 
for  any  finite  m,  but  fails  to  do  so  if  the  limit  m -*  ® is  taken  first. 
Conclusion:  square  of  wavefunction  amplitude  near  energy  minimum  at 

surface  vanishes,  due  to  b.c.  This  effect  disappears  deep  within  the 
bulk. 

The  case  (b)  is  more  complex.  Defining  6 s %(k^  - kQ)a  and  pro- 
ceeding to  the  limit  6 -»  0,  we  obtain: 

(b)  U -*  2N  ^ cos  k a sin  k am  sin  6m  (6  -*  0) 

z m m 

in  which  k a = cos  ^(-K/4L)  (cf.  Figure  b) . U.  and  U_  cannot  be  dis- 
ro  L z. 

tinguished  at  any  finite  value  'of  the  index  m (although  presumably  one 

of  the  two  solutions  has  an  extra  node  at  m = %N  as  compared  with  the 

z 

other,  to  ensure  orthogonality).  Note  that,  once  again,  as  we  approach 

2 2 

the  energy  minimum,  |u|  o,  now  proportional  to  6 .at  finite  m. 
Similar  results  obtained  at  the  minimum  of  (c) , and  at  the  maxima  in 
all  cases.  We  conclude  that  the  vanishing  of  the  square  of  the  wave- 
function  amplitudes  near  the  band  extrema  at  the  surface  of  a solid  is 
a general  feature  of  this  band  structure  (first-  and  second -nearest 
neighbor  overlap) . 

A Con  lecture: 

We  believe  that  the  above  conclusion  is  a general  feature  of  ar- 
bitrary band  structures,  required  because  of  some  as -yet -undiscovered 
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theorem.  A heuristic  proof  may  be  given  as  follows:  an  arbitrary  band 
structure  having  a minimum  at  k = 0 may  be  approximated  by  case  (a) 
near  the  minimum,  and  the  analysis  above  is  then  applicable.  An  arbi- 
trary band  structure  with  a minimum  at  k 4 0 may  be  approximated  by 
(b)  near  the  minimum,  and  the  analysis  appropriate  to  that  case  then 
applies.  In  both  instances,  the  square  of  the  amplitudes  |u|  -*  0 pro- 

protional  to  (E-€  . ) . Computations  carried  out  on  several  different 
min 

examples  of  band  structures  have  confirmed  these  heuristic  arguments. 
2D-Brillouin  Zone : 

For  each  (k^,k^)  there  is  one  k^  at  which  the  energy  is  a minimum 
and  a second  at  which  it  is  a maximum.  Denoting  these  by 
we  have  for  the  points  k]|  in  the  nth  band,  the  functions: 


f (k||,k  (kii)  . ) s € . (kji ) and  € (k,k  (kii)  ) = 6 (kii) 
nv  ,j  z II  ram  min  j|'  - n z U max  maxv  || 


defined  on  a 2D-Brillouin  Zone  spanning  all  allowed  values  of  kii 


(k  ,k  ) . This  is,  the  largest  possible  cross-section  of  the  3D-Bril- 
x y 

louin  Zone,  as  intersected  by  a plane  parallel  to  the  physical  surface. 


In  the  accompanying  paper,  Appendix  to  this  Chapter,  the  question  of 
the  stability  of  the  band  structure  against  the  formation  of  bound 
states  (surface  energy  bands)  is  treated.  It  is  seen  that  under  cer- 
tain conditions,  related  to  the  vanishing  of  a particular  component  of 

the  bulk  inverse-effective-mass  tensor,  a =0,  surface  states  are 

zz 

formed  even  for  infinitesimal  surface  perturbations.  The  Brillouln 
Zone  of  surface  states  may  cover  only  a small  portion  of  the  maximal 


2D-Brillouin  Zone.  For  large  surface  perturbations,  a bound  state  will 
form  below  each  energy  minimum  or  f°r  repulsive  perturbations, 


I 

I 


above  each  an^  the  Brillouin  Zone  of  bound  states  will  coin- 

'cide  with  the  maximal  size  of  a 2D-Brillouin  Zone.  The  reader  is  re- 
ferred to  the  Appendix  of  this  Chapter  for  examples  of  the  calculational 
method . 

Ad -A toms  on  Surface; 

We  consider  the  effects  of  an  atomic  perturbation  on  a surface  in 
the  presence  of  surface  energy  bands.  A perturbing  potential  always 
has  a bound  state  in  1 or  2 dimensions;  however,  the  bound  state  here 
may  coincide  with  the  continuum  of  bulk  states,  and  is  broadened  into  a 
resonance.  This  leads  to  the  possibility  of  the  following  possible  en- 
ergy levels: 


The  bulkd.o.s.  is  shaded,  surfa*ce  states  cross-hatched,  atomic 
level  dashed  line,  (a^  is  broadened  atomic  level, (b\ is  bound  state  below 
surface  band  (broadened  by  resonance  with  bulk  band)  and  (c) is  bound 
state  above  both  bands,  an  unbroodened  sharp  level. 
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APPENDIX  TO  CHAPTER  III 

ELECTRONIC  INSTABILITY  OF  SURFACES  OF  SOLIDS* 

Daniel  C.  Mattis 

Belfer  Graduate  School  of  Science 
Yeshiva  University 
New  York,  N-Y.  10033 


Abstract: 

We  study  the  effects  of  a surface  perturbation  on  a semi-infinite 
solid.  The  3D  energy  band  structure  is  found  to  determine  whether  or 
not  there  is  an  intrinsic  instability  against  the  formation  of  surface 
bands.  A criterion,  involving  one  relevant  component  of  the  inverse- 

m 

effective  mass  tensor,  is  derived. 


*This  research  supported  by  a grant  N00014-76C-0690  of  the  Office  of 
Naval  Research.  's 
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Experimental  and  theoretical  analyses  have  by  now  established 

3 

that  the  surfaces  of  solids  often  differ  in  crystallographic  class  as 

4 

well  as  in  electronic  properties  from  the  bulk.  The  fact  that  not  all 
surfaces  have  been  found  to  be  reconstructed  led  me  to  seek  a criterion, 
and  ultimately  to  pose  a simpler  problem:  under  what  circumstances  will 
an  infinitesimal  perturbation  produce  a band  of  surface  states?  We 
shall  denote  this  an  intrinsic  instability,  for  while  any  arbitrary 
material  may  or  may  not  have  surface  states,  those  with  intrinsic  instabil- 
ilities  always  must?  In  the  present  work,  we  shall  prove  that  a surface 
lying  in  the  x-y  plane  is  intrinsically  unstable  if  and  only  if  a compo- 
nent of  the  bulk  inverse-effective  mass  tensor  or  vanishes  at  an  appro- 
priate point  or  set  of  points  in  the  Brillouin  Zone.  It  is  thus  the  bulk 
band  structure  and  the  surface  orientation  that  primarily  predetermines 
the  intrinsic  surface  instability  --  the  nature  of  the  surface  perturbation 
potential  is  of  secondary  importance. 

Our  prototype  material  consists  of  atomic  planes  at  z = na,  with 
n = 1,2,...,N  the  plane  index.  It  is  terminated  by  the  two  surfaces  at 

n = 0,  N +1,  with  the  x and  y coordinates  continued  periodically  at  N +1 
z X 

and  Ny+1  respectively.  Proceeding  to  the  limit  and  all  -*  °°,  we 

concentrate  on  the  one  surface  at  n = 0.  Consider  the  "unperturbed" 
Hamiltonian  Hq,  having  the  following  matrix  structure  within  a single 


Bloch  band; 


I V C,<')  H„  e M-  €.  K 

k'  k k k,k' 


(1) 


in  which  *s  the  appropriate  linear  combination  of  Bloch  functions  Qf^, 
which  satisfies  the  boundary  condition  = 0 for  s < 0.  In  the  case 
of  a simple  band  structure  having  (at  fixed  k ,k  ) a minimum  att  k a ■ 


x y 


or  tt,  this  linear  combination  is  merely: 


and  the  first  layer  amplitudes  are  Cj<(x»y»a)  = X^(x ,y)f (kza)  with 


f(k  a)  = sin  k a, 
z z 


In  cases  when  the  band  structure  €(k  tk  ,k  ) has  a minimum  at 

x y z 

kza  ^ 0 or  it  (with  k|j  = (k^,k^)  f ixed) , the  amplitude  near  the  surface 
may  be  obtained  from  a solution  of  the  SchrHdinger-like  equation: 


€(kii,ij£)  f(z)  = CCkiijk^  f(z),  for  z > 0, 


subject  to  the  boundary  conditions  f(z  < 0)  =0.  We  have  established 

that  for  k near  the  energy  minimum  fee  sin(k  -k  . )z,  and  obtained  a 
z z z min 

similar  result,  f CCsin  (k  -k  )z,  near  the  energy  maximum.  These  are 

Z ZUlclA 

basically  the  only  properties  we  shall  require. 

We  next  assume  that,  either  as  a result  of  a nonzero  change  in 

Madelung  potential"*  at  small  values  of  z,  or  because  of  a small  displace 

ment  of  the  surface  atoms  from  their  ideal  positions,  or  a result  of  any 

other  physical  requirement,  a perturbation  Hamiltonian  H exists  near 

s 

the  surface.  Explicitly  factoring  out  the  amplitudes  f(kza)  of  the 
scattered  waves,  as  dictated  by  geometric  considerations,  we  may  without 
loss  of  generality  write  the  matrix  structure  of  Hg  in  the  form: 


■f  v £<?)  v?)  y'>=  < wy  6kii,k'ii  f<k>>  <4> 


diagonal  in  kj|  because  of  translational  invariance  in  x and  y.  If  the 
perturbation  is  restricted  to  the  first  surface  layer  (z=a)  the  coupling 
constant  is  a function  only  of  k jj , and  is  denoted  g(kj|).  If  the  pertur- 
bation affects  several  surface  layers  the  coupling  constant  will  depend 


also  somewhat  on  k^  although  for  simplicity  In  the  ensuing  calculations 

we  approximate  it  by  itr.  value  at  ^zmjn^or  attractive  perturbations)  or 

k (repulsive)  so  that  g.  . , is  then  just  a function  of  kii,  again 
zmax  k,k  || 

denoted  g(kj|).  In  fact  our  criterion  for  inherent  instability  does  not 
depend  on  the  value  of  g,  merely  on  its  sign  and  on  the  fact  that  at 
^zmin/max  t*ie  couP^^n6  constant  is  neither  identically  zero  nor  infinite 
Ultimately,  if  our  procedure  is  found  to  lead  to  interesting  results,  it 

g 

may  be  rendered  quantitative  by  obtaining  g^  self-cons  is tently  . 

We  first  study  the  bound  states  of  the  joint  Hamiltonian  H + H , 

J J os 

l 4 

which  if  they  exist,  constitute  the  much  discussed  ’ "surface  energy 
bands".  They  must  be  of  the  form 


♦kiu-v’  *'<kl-V  'kll.k. 


(5) 


in  which  we  use  "s"  as  a subscript  to  distinguish  these  surface  modes 

from  the  continuum  of  bulk  states.  Schrildinger 's  equation  (HQ+Hg) 

«*  E y yields  the  following  equation  for  the  coefficients  F: 
s s 


(€  -Es)P(k||,kz)  + i £ - o (6) 


k' 

z 


In  order  for  this  equation  to  have  a nontrivial  (F  4 0)  solution,  the 
following  secular  equation  must  be  satisfied: 


<k||>  - 


(7) 


in  which  ^ is  a special  case  of  a more  general  function: 


3 


t ; 


M’ 

l ; i 
' 

E I 


on  the  surface  orientation  which  defines  the  z-axls  relative  to  the 
invariant  crystal  axes.  It  is  convenient  to  write  the  bound  state  energy 
in  the  following  analogous  form: 


€k||>kmin+  °U(k||-'W„>  (1-cosh  X) 


Oil) 


Now,  Eq.  (7)  simplifies  to: 

- ± a <k.,,k  , ) = i f d0  e‘IXl  H2) 

g zz  ||  zmin  tt  cosh  X + cos  0 

o 

Thus  the  r.h.s.  of  Eq . (12)  fails  to  diverge  when  \ - 0 (in  fact, 
approaches  a maximum  limiting  value  of  l) . There  may  nevertheless  be 
an  incipient  instability  against  the  formation  of  surface 'states  at 
arbitrarily  small  |g  | provided  there  exists  a locus  in  2D  k||-space  for  which 
azz  ■ 0*  Examples  of  this  include  the  elementary  tight-binding  s-bands 
* -K(cos  k a cos  k a cos  k a)  for  b.c.c.  and  -K(cos  k a cos  k a 

k x y z * y 

+ cos  k a cos  k a + cos  k a cos  k a)  for  f.c.c.  crystal  structure,  both 
having  a minimum  at  k^a  = 0 or  tt.  For  each  there  are  rectilinear  segments, 
along  which  (k||,0)  = 0.  Therefore,  (12)  will  have  solutions  along  a 
neighborhood  of  these  lines,  over  a surface,  the  area  and  geometry  of  which 
will  be  functions  of  the  magnitude  of  g.  Conversely,  for  the  s.c.  tight- 

m 

binding  band  structure  -K(cos  k a + cos  k a 4-  cos  k a)  for  which  the  minimum 

x y z 

is  also  at  k = 0,  there  exists  no  curve  over  which  or  ■ 0,  and  therefore 

Z zZ 

no  region  of  ky-space  inherently  unstable  against  the  formation  of  surface 
states. 

The  analysis  for  g > 0 is  similar.  The  Bloch  energy  is  expanded  about 
its  maximum,  the  effective  mass  parameter  is  now  non-positive,  and. the  bound 
state  lies  above  the  continuum.  With  the  aid  of 
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£ (E 
n,m 


, f(k'na)  f(k'ma) 

;X||)  ■ h:1  L I — 5 5 — I 


k‘  E - «kll,k"l0+ 


n/j 


* fo  ^ t(-K  “>  f<k>>  + 1"«  <E-\||,k.)  I «) 


Normalization  of  the  wave  function  requires: 


- £<k.a) 
F(k) 


E - € 
s 


c ~a  Sl.l^lP] 


3E 


(9) 


As  Eg  approaches  €(k j|  kzm/n)  from  below  (for  attractive  potentials, 


or  as  E approaches  €(kn,k  ) for  repulsive)  the  denominators  become 
s 1 1 z max 


arbitrarily  small.  However,  the  numerators  also  vanish  (cf.  discussion 


following  Eq.  (3))  and  thus  the  integrals  remain  finite  unless  — due 


to  the  peculiarities  of  the  band  structure  — the  denominator  vanishes 


as  (k  ~k  ) . The  conditions  for  this  to  occur  are  now  examined.  For 
z zm 


g < 0,  we  start  by  expanding  the  Bloch  energy  about  its  minimum: 


^ll-k ^U.W 4 ^.<k||-W<Vkz»i/  »2-  °«vw4> 


or,  more  compactly,  with  0 = (k  -k  . )a, 

z zmin 


’k-"  V’Wn*  * z(k||’kzmin)a-COS  9)  + ° 


(10) 


In  which  orzzOc||»k2mtn)  H (®  )yz  s a 2d2€_4  /Sk2  is  the  zz-component  of  the 

k 


zz  ' ||'  zmin  'zz 

inverse-effective-mass  tensor.  Because  it  is  evaluated  at  a minimum. 


o ^ 0 necessarily;  its  magnitude  depends  both  on  the  band  structure  and 

ZZ 


I 
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~ V'W'  “ S(k||-kama«>(l  + c°°  0>  + 0 <**> 


"z  \n,k  4-  a (kn,k  ) (1-cosh  X) , 
||  zmax  zz  ||  zmax 


Eq.  (7)  once  again  reduces  to  (12),  after  the  substitution  of  Of  (kii.k  ) 

zz  ||  zmax 

for  °zz  (^J|»^zmin^  *-n  t^le  Because  of  the  change  in  sign  of  the 

effective-mass  parameter,  the  bound-state  solutions  now  exist  only  for 

g > 0,  but  aside  from  this,  the  discussion  given  after  Eq . (12)  applies 

for  this  case  also.  It  should  be  noted  that  the  loci  of  a =0  are 

zz 

not  necessarily  the  same  for  k as  for  k . , although  in  the  two 

zmax  zram 

examples  given  (b.c.c.  and  f.c.c.)  they  do  coincide  because  of  symmetry. 

The  s.c.  tight-binding  band  structure  provides  an  example  of  the 

importance  of  surface  orientation.  Inherently  stable  against  the 

formation  of  surface  bands  when  the  surface  is  a (100)  plane,  it  demonstrates 

an  inherent  instability  along  a (Oil)  direction:  in  the  new  coordinate 

system  a=a  =a2  ,a=a,  and  the  Bloch  energy  takes  the  form 

y z x 


€ = -K  (cos  k a + 2cos  k a • cos  k a ) 

V x y y z z 


For  \maz~  0,tt  the  equation  a ^ = 0 has  a solution  along  the  straight-line 
segments  k^  = arbitrary,  k^  = + rr/2a^,  indicating  an  instability  against 
the  formation  of  surface  states  in  the  neighborhood  of  these  lines. 


A surface  perturbation  scatters  bulk  states,  and  modifies  the  surface 

density  of  states  g affecting  optical  absorption  at  the  surface,  etc., 

even  when  there  are  no  states  bound  to  the  surface.  The  contributions  of 

bound-  and  scattered-  states  are  combined  in  the  calculation  of  the 

2 

"perturbed"  local  density  of  states  function  on  the  n-th  plane.  After  some 
algebra,  one  obtains: 


LI  l 


* J 


! 


• 1 1 ■■ — — — — - - 


»n(«i)  - oy./1  l!  <^>n<vk|l»  <:ah,i<Vkll)/a!.)’l*(,,rE.) 


II 

t > w1  c IV  n(k)|2  6(oh €J 


(14) 


where  E is  the  surface-band  energy,  and 


(-  _f(%»a)Cq-g(k||)S1-t(^r;k||)>g(ic||)f(k8a)Sl_ii(6|?;l.||) 

1 • E(k||)  El,l(ej:k||)  (15) 


In  the  bulk  limit  n -*  00  we  drop  the  exponentially  decaying  and/or  rapidly 


-1 


oscillating  terms  to  obtain  the  limiting  behavior  p (id)  = (N  N N ) £ 

00  V U 7 “* 


x y z 


6 ( 03  — (:_) , which  is  exactly  the  usual  bulk  value.  Thus,  the  existence 
k 

of  the  surface  is  "forgotten"  deep  within  the  bulk.  For  small  n 


however,  the  local  dens lty-of-states  given  by  (14)  is  sensitive  to  the 


sign  and  strength  of  the  perturbing  potential,  and  to  the  existence  of 
10 


bound  states. 

The  modified  charge  distribution  is  conducive  to  the  phenomenon  of 

3 8 

"surface  reconstruction"  * . We  have  seen  that  the  surface  states  (if  any) 
do  not  necessarily  occupy  the  full  2D  surface  Brillouin  Zone,  but  just  a 


small  area  near  the  locus  of  the  instability  Of  = 0.  The  narrowness  of 

z z 


this  area  in  k-space  suggests  a large  resultant  cell  size  in  real  space, 
as  is  indeed  experimentally  observed  in  some  reconstructed  surfaces.  The 
predicted  absence  of  inherent  surface-state  instability  in  other  cases 
(e.g.  the  model  s.c.  with  (100)  surface)  suggests  the  absence  of  surface 


reconstruction  there.  In  all  cases, an  analysis  of  a (k n , k ) to  determine 

zz  zm 


where  (or  if)  O’  = 0 might  usefully  precede  any  full-blown  study  of  a sur- 
ZZ  • 


face,  to  pin-point  the  a-priori  instabilities. 

I thank  Drs.  Rafael  Pena  and  Barry  Simon  for  helpful  collaboration. 
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We  have  completely  solved  two  examples:  € = ( — ;r  4-  k . ) , with 

dz2  zmin  ’ 

the  helpful  collaboration  of  Prof.  Barry  Simon,  and 

€ = -A  cos  a + B cos  2a  . In  both  instances  f(z)  has  the 
behavior  indicated  in  the  text,  and  we  have  compelling  reason  to 
believe  that  its  vanishing  near  k . ( or  k ) is  a very 
general  property;  a full  account  will  be  published  elsewhere. 


Mathematically,  the  problem  of  obtaining  bound-  and  scattering-states 

of  H + II  remains  exactly  soluble  for  any  g ,,  which,  while  an 
os  k , k 

arbitrary  function  of  ky,  is  of  the 

separable  form  0(kz)0(kz ')  in  k^ , or  is  a finite  sum  of  each  terms. 

Instead  of  (7)  we  obtain  a secular  determinant  containing  the  quantities 

E . The  condition  for  the  existence  of  a bound  state  when  |c  I-*  0 
n,m  1 

remains  precisely  a = 0 in  this  case.  In  the  most  general  case, 

zz 

g , , an  arbitrary  function  of  k and  k',  we  once  again  recover  the  same 
Kj  K Z Z 

criterion  using  a variational  solution  for  the  bound  state.  This 


requires  only  that  g^  be  suitably  bounded  — neither  zero  nor  infinite 
at  k 

zm 

8.  As,  for  example,  in  the  work  of  M.  Schldter,  J.R.  Chelikowsky,  S.G.  Louie 
and  M.L.  Cohen,  Phys.  Rev.  Letters  34,  1385  (1975).  In  our  method, 

the  calculation  of  a self-consistent  surface  perturbing  potential  Hg 
would  involve  the  perturbed  dens ity-of-states  Eq.  (14)  and  Poisson's 
equation. 

9.  as  observed  in  photoelec tronic  spectroscopy:  see  D.E.  Eastman  and 
M.I.  Nathan,  Physics  Today  2jJ,  44  (1975)  and  R.L.  Park,  ibid . . p.  52. 
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10.  The  reader  should  note  that  the  continuum  of  surface  states  in 

general  overlaps  the  bulk  band  continuum  of  Bloch  energies.  Thus 
the  surface  bands,  if  any,  may  lie  entirely  within  the  Bloch-state 
band  (this  is  the  case  of  an  intrinsically  unstable  surface  under 
the  influence  of  a very  weak  Surface  perturbation),  or  may  only 
partly  overlap  it  (somewhat  stronger  perturbing  potential)  or 
may  indeed  lie  entirely  outside  the  bulk  band  (as  happens  when  the 
perturbing  potential  is  large, regardless  of  whether  or  not  the 
surface  is  intrinsically  unstable  in  the  sense  of  the  present  work.) 


IV.  TRANSPORT  IN  A RANDOM  MEDIUM 


Impurity -band  conduction  and  electron  energy  bands  in  random  al- 
loys may  be  studied  by  similar  techniques.  The  best  single  method  of 
the  past  decade  has  been  the  CPA,  which  suffers  from  a basic  weakness 
inability  to  take  short  range  order  (SRO)  into  account.  A relatively 
simple  method,  yet  one  that  can  account  for  SRO,  and  LRO  if  necessary 
is  outined  in  the  Appendix  to  this  Chapter.  Extension  of  this  to  the 
Hubbard  model  (2 -body  forces)  is  actively  underway,  the  SRO  being  de- 
termined self-consistently . 
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APPENDIX  TO  CHAPTER  IV 


ELECTRON  STATES  IN  RANDOM  ALLOYS  WITH  SHORT-RANGE  ORDER 


* it 

Paul  Bloom  and  Daniel  Mattis 


Belfer  Graduate  School  of  Science 


Yeshiva  University, 
New  York,  N-Y.  10033 


Abstract: 

We  present  an  accurate  and  economical  iterative  method  of 
calculating  the  energy  levels  of  a disordered  or  partly  ordered 
random  alloy.  Results  presented  for  1-D  and  3-D  simple-cubic, 
lattices  compare  favorably  with  exact  calculations.  We  also 
present  the  systematic  effects  of  partial  short-range  order  in 
3-D.  A theory  of  the  one-particle  propagators  is  presented, 
and  the  theory  of  electrical  conductivity  is  developed  in  the 
context  of  our  new  method.  Our  formulas  satisfy  the  exact 
conservation  laws. 


* Supported  by  AFOSR  Grant  it 73-2430B- 


it  This  research  is  supported  by  a grant  from  the  Office  of  Naval 
Research  #N0001A -76-C-0690 
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I Introduction 
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The  study  of  electronic  and  vibrational  spectra  of  disordered 

alloys  is  currently  one  of  the  principal  concerns  of  solid  state 

physics  , stimulated  by  the  outstanding  successes  of  the  coherent 

2 3 

potential  approximation  ’ (CI'A),  now  ending  its  first  decade. 

Although  efforts  to  improve  our  understanding  beyond  the  CPA  have 

not  all  met  with  the  same  good  fortune,  there  have  been  recent 

4 5 

exceptions.  Cluster  methods  ’ have  been  devised  which  are  accurate 
enough  to  reproduce  the  "peaky"  structure  of  the  density  of-states 
P ( a>> , which  they  sometimes  do  (notably  in  one  dimension^)  with 
startling  fidelity.  We  have  been  working  along  such  a cluster-type 
approach,  and  have  found  an  extremely  simple  method  translating 
directly  into  a computer  algorithm.  While  unsuited  to  the  theoretical 
studv  of  Lifshitz  tails,  our  method  has  permitted  us  to  reproduce 
many  of  the  other  known  results  over  the  theoretically  permitted 
range  of  energy^,  even  near  the  energy  maxima  and  minima,  and 
additionally,  permitted  to  study  of  the  effect  of  short-range  order 
(SRO)  . Along  with  Lifshitz,  we  envisage  tails  in  p(iu)  at  the  energy 
maxima  and  minima  as  arising  from  accidental  correlations  in  increas- 
ingly large  clusters,  of  a size  that  for  practical  reasons  we  are  not 

at  present  capable  of  handling;  however,  the  simplicity  of  the  present 

8 

method  may  suggest  a natural  extension  to  cover  this  . 

The  basic  outline  of  our  paper  is  as  follows:  In  Section  II 
we  present  a method  for  the  calculation  of  the  single-body  Green 
function  in  the  presence  of  an  arbitrary  number  of  impurities.  We 
then  discuss  how  our  procedure  can  be  implemented  by  the  use  of  a 
convergence  factor,  £.  Section  III  is  devoted  to  an  analysis  of  the 
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meaning  and  uses  of  the  complex  self-energy  Z within  the  context  of 

a disordered  medium.  Results  from  our  method  are  presented  in 

Section  IV,  including  the  effects  of  SRO.  Beyond  this  in  Section  V we  make 

further  approximations  that  allow  us  to  determine  _^( oj)  . 

k k 

Section  VE is  concerned  with  the  development  of  a transport  theory 

9 

compatible  with  G,  along  the  lines  of  Baym  and  Kadanoff. 

II  Cluster  Green  Function 

Let  the  Hamiltonian  for  the  electrons  within  a single  tight 
binding  band  in  a hypercubic  lattice  in  D dimensions  be: 


H = £ T.  . 1 i>< j 1 + Z V |i><i|  = T + V 
i,j  i 


with  T.  . ■=  (2D)  * for  i,j  nearest-neighbors  and  zero  otherwise, 

Ji)  the  Wannier  state  at  the  lattice  point  and  V^^  the  potential 
which  takes  on  one  of  two  values  depending  whether  atom  A or  B 
occupies  the  i site.  We  construct  the  resolvent  operator  G(z)  and 
its  various  matrix  elements: 


G(z)  = (z-H)' 


[z  - (T+E)  - (V 


in  which  we  reference  the  operators  to  a complex  "optical  potential" 
E(z). merely  as  a device  to  enhance  the  convergence  of  subsequent 
expansions,  with  z = frequency  w, extended  to  the  complex  plane. 
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For  those  readers  familiar  with  the  CPA,  it  is  important  to 
note  that  our  new  departure  consists  principally  in  dissociating  the 
complex  self-energy  parameter  £(u>)  from  the  site-diagonal  averaged 
Green  function  Gnn(iu)  . Whereas  in  CPA,  knowledge  of  the  one  implies 
the  other,  via  the  relationship: 

(*)  = <n| [ 2 - (T  + £CPA(z»:f1|n>  > (3) 

CPA 


our  experience  indicates  that  it  is  better  to  treat  2(a))  merely  as  a converge 

parameter,  one  to  be  chosen  as  an  ad-hoc  aid  in  the  calculations 

rather  than  by  tedious  and  unnecessary  self-consistency  conditions.  As  by 

r.n . (21  the  exact  G ( z)  are  all  indpnendpnr  n'f  z) , in  apv  accurate  apnrex 
^ v * nm  * * 

mation  to  G (z)  we  have  latitude  in  our  choice  of  £(z),  as  discussed  below 
— — — nm  \ * 

and  we  pick  the  simplest  possible  £(z)  for  which  our  calculated  G is 
approximately  stationary. 


We  next  define  a modified  resolvent  operator,  G^\  appropriate 
to  the  case  in  which  one  sets  = 0,  where  we  define  Vi  = <i|(V-£).|i) 
and,  indicating  the  elimination  of  the  localized  fluctuation  potential 
at  this  site  by  ( )^,  we  have: 


G(l)(z)  = [z  - (T  + £)  - (V  - E)  'J'1  . 


The  full  resolvent  (2)  can  be  expressed  in  terms  of  the  modification 
in  (A)  by  the  use  of  the  operator  identity  (A-B)'1  = A-1  + A*1  B(A-B)"1; 

G(z)  « G(i)(z)  + G(i)(z)  (V  - £).  G(z)  . (5) 

55 

*■  ihAhMS,  , |ffl  j 1 ...  — - I.  


- 


1 £ 


\ K 

i 


Because  the  perturbation  is  diagonal  in  the  Wannier  representation,  the 
matrix  elements  are  easily  found: 

m Cni><z>  Vi  Gi^<2> 

G«»(z>  ■ <2>  + — (77— r • 

1 - G)y(z)  V. 

u i 


For  the  calculation  of  the  density  of  states  function  p(u>)H  — ImGnn(iu+ie  ) 
only  the  configurationally-averaged  Gnn(io+ie)  is  required. 

For  the  one-particle  propagators  g the  averaged  Fourier  transforms 



of  all  (u>  + ie)  are  needed.  Eq.  (6)  is  now  iterated.  Define  G^’^Cz) 
to  be  the  modified  resolvent  operators  with  the  fluctuation-potentials  V 
at  sites  i and  j removed.  By  a repetition  of  the  above,  we  have: 


G^<z) 


G^1* J)(z)  + 

nm 


,G<r’j)  V.  G^1’^ 

nj  j jm 

1 - g\x*^  V. 
jj  J 


The  matrix  elements  G decay  exponentially  with  distance  R : 

nm  nm 

thus  the  expansion  (6),  (7)  is  in  a symbolic  "parameter"  y defined  as 
Gn/’’’^Vj»  which  is  "small"  for  small  and  "exponentially 
small"  at  large  V ^ v The  process  (6),  (7)  is  to  be  repeated  any 

m 

numer  of  times,  until  fhe  largest  practical  cluster  size  is 
10 

achieved . Termination,  by  truncation f of  the  series  consists  of 
approximating  the  most  distant  G's,  i.e.  those  with  the  largest 
number  of  superscripts,  by  their  value  in  the  average  optical 


potential. 


Thus,  if  we  stop  at  (7), 


the  approximation 


consists  in  replacing  G^,’^(z)  by  <n|[z  “ (T  + S)]-1|m).  The  configurational 

averages  over  all  the  explicitly  retained  are  then  performed,  and  all  G's 
obtained. 
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Ill  Choice  of  £ 

We  now  come  to  our  principal  point  of  departure  from  other 
methods  - our  choice  of  £.  Our  results  would  depend  crucially  upon 
£ except  for  the  following  observations.  Since  the  behavior  of  the 
local  cluster  is  the  dominant  characteristic  of  disordered  systems, 
we  expect  results  insensitive  to  the  particular  choice  of  £ if  the 
cluster  size  is  sufficiently  large. 

We  require  a simple  functional  form  for  £ that  allows  for  states 
out  to  the  bands  limits.  This  excludes  the  use  of  £ . which  is  known 

Ll  A 

to  produce  bands  that  are  always  too  narrow.  We  restrict  the  range 
of  possible  £'s  by  requiring  that  it  obey  dispersion  relations,  insur- 
ing that  our  approximate  G is  analytic.  Furthermore,  a functional 
form  is  desired  in  which  G is  accurate  in  both  the  weak  as  well  as 
strong  scattering  regimes.  Because  of  the  local  nature  of  highly 
disordered  systems,  our  choice  becomes  more  critical  for  small  potential 
differences  where  effects  are  more  extended.  Our  input  is  the  Im  £ 
which  we  take  as  one  or  more  step  functions , non-zero  only  within  the 
theoretical  band  limits.  Re  £ is  then  determined  from  the  following 
dispersion  relation 


E(z) 


t<vL  + i_  r 

2 +2TTi  J 


dx 


££x+ie) 

x-z 


(8) 


This  is  sufficient  to  make  our  approximation  to  G(z)  satisfy  causality. 

00 

The  density-of-states  sura  rule,  J duP(ili)  ° l is  itself  a beneficial 

-00 

consequence  of  the  analyticity  of  our  approximate  G(z)  and  its  result- 
ing l/z  dependence  in  the  asymptotic  limit  as  we  discuss  elsewhere. 

We  verified  that  in  all  cases  studied,  the  sum  rule  on  p(u>)  was  satisfied 

« 

numerically. 
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If  we  choose  the  constant,  Im  £,  to  be  of  magnitude  of  Im  £ 

t*  A 

then  in  the  weak  scattering  and  low  concentration  regimes  G(z)  will 
be  quite  similar  to  GcpA(z)  so  that  accurate  results  can  be  expected 
in  all  regimes. 

Before  we  proceed,  the  way  in  which  we  use  Im  £ must  be  more 

CPA 

clearly  outlined.  In  the  accompanying  Figure  l we  display  the  two 

12 

basic  behavior  patterns  of  |lm  £cpA  | as  observed  by  VKE . It  should  be 
noted  that,  here  also,  Re£  arx!  Im£  are  related  by  the  Eq.  (8).  £ has 

to  describe  everything  in  the  CFA;  it  determines  band  gaps,  peaks  in 
the  density-of-states , and  the  general  overall  scale.  Most  of  these 
results  (e.g.  band  gaps  and  complicated  structure)  are  better  obtained 

by  our  detailed  calculations  of  the  correlated  scattering.  We  hypothesize 
that  the  most  useful  information  from  CPA  is  contained  in  the  general 
overall  magnitude  of  Im  ^CpA*  Operationally,  in  Figure  la,  we  would 
. ignore  values  of  |lm  £cpA  | from  the  region  of  its  maximum  as  well  as 
the  extreme ties  of  the  band.  In  the  former  range  of  energies,  we  expect 
exceptional  scattering  because  it  is  easiest  for  these  states  to  make 
transitions  due  to  band  overlap,  whereas  at  the  band  edges  the  spectrum 
will  be  least  disturbed,  according  to  the  same  considerations.  Any 
value  from  the  shaded  region  is  then  acceptable.  In  terms  of  particle 
lifetimes,  we  will  obtain  the  large  and  small  transition  rates  because 
we  almost  solve  the  eigenvalue  problem  exactly  for  each  configuration 
and  this  is  clearly  equivalent  to  a perturbation  approach.  As  for 
concentration  dependence  in  we  will  obtain  correct  behavior  simply 

because  we  weigh  each  configuration  by  its  appropriate  probability.  Thus 
we  are  able  to  include  both  the  dynamical  and  statistical  aspects  of  the 
problem. 
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In  case  b,  the  same  analysis  leads  us  to  ignore  the  very 


••t 
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f 
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large  values  of  jlm  E | in  both  subbands.  Here  though,  the 
magnitudes  are  considerably  different  leading  us  to  suspect  that  two 
different  constants  are  needed.  Further  details  of  this  case  will 
be  elucidated  in  the  following  examples. 

IV  Analysis  of  Results 

We  first  consider  the  canonical  l-D  tight  binding  binary  alloy  for 
3 different  scattering  strengths  at  a 50-50  concentration. Figure  2 
compares  the  results  of  1,3,  and  5 cluster  calculations  for  p(ir)  when 

= + .5  with  exact  results.  We  see  in  this  example  the  development  of  the  peaky 
structure  associated  with  special  clusters  of  atoms  as  our  cluster  size 
increases.  Proceeding  to  a larger  scattering  strength  (V^  = +1.0), 
we  expect  that  the  local  configurations  will  play  a more  prominent 
role  because  of  increased  wave  function  localization.  As  shown  in 
Figure  3 we  successfully  reproduce  most  of  the  structural  details  of 
p(u>)  for  a 5 cluster.  To  check  the  degree  of  insensitivity  in  our 
5 cluster  model  we  varied  |lm  £|  within  the  limits  given  by  Im  I 

CrA 

and  found  little  change  in  the  overall  pattern  as  shown  in  Figure  4 
This  indicates,  numerically  that  the  resulting  G is  stationary  and 

m 

that  E is  optimum. 

The  acattering  strengths  arc  now  increased  to  = + 2.0, 
providing  a critical  evaluation  of  the  methods  capabilities  (larger 
scattering  strengths  are  in  a sense  too  easy  because  wave  function 
localization  makes  a cluster  calculation  more  plausible).  Using  the 
exact  scattering  off  all  configurations  of  5 atoms,  the  highly  discrete 
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spectrum  is  well  reproduced,  as  seen  in  Figure  5a.  Increasing  the 
cluster  size  to  7 atoms,  keeping  the  convergence  factor  the  same  as  in 
5A,  improves  our  agreement  with  the  exact  results  as  shown  in  Figure  5b. 
In  Figure  6 we  display  the  results  of  again  varying  |lm  £|  within  the 
limits  dictated  by  Im  £ ; the  major  details  are  again  seen  to  remain 

v*n 

stationary.  We  have  found  empirically  that  if  |lra  £|  is  too  small,  the 

resultant  density-of-states  is  too  "peaky"  and  as  such,  representative 
of  a molecular  cluster,  instead  of  the  solid  state.  If  |lm  £|  is  too 
large,  then  the  central  site  predominates,  as  is  correct  only  in  the 
extreme  "atomic"  limit  when  potential  fluctuations  greatly  exceed 
the  band-width.  One  can  see  this  from  Figure  6 since  the  sharper 

curve  is  associated  with  the  lowest  value  of  |lm  £|  and  vice-versa. 

In  three  dimensions  the  obvious  cluster  size  is  7 sites.  Figure  7 

compares  our  calculation  with  the  Monte-Carlo  type  numerical  results 
13 

of  Alben  ejt.  al_.  A constant  Im  £ gave  poor  results  in  this  case,  but 
the  CPA  calculations  immediately  showed  us  why  ! Im  £ was  more 

than  one  order-of-raagnitude  smaller  in  the  majority  subband  than  in 
the  minority  subband.  Consequently  we  changed  Im  £ to  the  step  function 
shown  in  the  figure,  varying  the  parameters  (magnitudes  of  the  steps) 
again  guided  by  CPA.  The  results  now  agreed  well  with  the  exact 
computations  and  were  insensitive  to  the  precise  value  of  our  parameters 
as  is  evidenced  by  Figures  in  which  a 3-step  function  was  used. 
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To  illustrate  entirely  new  applications,  consider  effects. 

of  SRO  on  this  same  alloy.  With  the  Cowley  SRO  parameter,  c 

A 

and  cfi  = 1 - cA  the  relative  concentrations  and  PAB  the  probability 
of  finding  atom  A at  a given  site  when  a B atom  occupies  a specified 
neighboring  rite,  „e  have  P^  =■  cA  + c,  n,  PM  - cB(i  - a),  ?AB'-  eA(l  - „)  , 

Snd  PEB  ■ CB  + CA°'-In  tef-  13.  o-o.  For  cA  » .1,  ct  can  vary  from 
-0.11  to  +1.0;  negative  a is  associated  with  enhanced  tendency  of 
A atoms  to  be  surrounded  by  B’s  (i.e.  "antif erromagnetism") , positive 
a indicates  enhancement  in  the  probability  of  either  species  being 
surrounded  by  atoms  of  its  own  kind  (i.e.  "ferromagnetism”).  Using 
the  same  convergence  parameters  as  in  our  calculation  at  = 0,  In  Pigure  9 
we  find  distinctive  features  in  the  minority  sub-band  density  of 
states  that  we  interpret  in  terms  of  minority-atom  clustering: 

the  single  peak  of  q = -.07  registers  cne  unlikelyhood  of  finding 
two  A atoms  as  nearest  neighbors,  and  the  double  peaks  of  a = .7 
represent  the  tendency  of  the  same  atoms  to  form  pairs,  triplets , etc. 
However,  due  to  the  sparseness  -of  A atoms,  triplets  and  higher-order 
clusters  are  statistically  insignificant  for  these  values  of  a . 

V Electron  Propagation 

So  far  we  have  developed  a 1 ethod  for  calculating  the  site-diagonal 
configuration  averaged  Green  function.  We  have  not  indicated,  how  we 
would  calculate  the  non  site-diagonal  propagators.  One  alternative  is 
to  develop  a cluster  method  for  the  latter,  similar  to  the  method  we 
used  for  the  former.  Another,  simpler  though  less  accurate,  alternative 
will  be  employed.  We  first  define  a new  self-energy  A*(z)  by  the  equation 
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(9) 


V2> ■ » 5 


z - A*(z) 


This  relationship  is  numerically  inverted  to  obtain  A*  as  a function  of 


the  exact  or  numerically  calculated  G^  . From  the  calculation  of 
A*(tu)  we  obtain  A*(z)  in  the  entire  complex  plane,  as: 


A*(z)  = <V)  + ^ ^ dm  Im  A*(m+iR). 

0)-Z 


(10) 


If  further  values  of  G (z)  were  calculated  numerically .then 

AJ 

we  would  determine  A*(k,z)  from 


G 


(z) 


and 


. — * — * v 
ik.(R  -R  ) 
e . H m 

z - A*(£,z)  - e 

k 


(ID 


1 . 

z - A*  (k , z ) - e_^ 
k 


The  analysis  is  facilitated  by  going  over  to  a localized  representation 

in  which  we  would  specify  the  number  of  elements  A*  (z)  that 

£m 

we  have  determined  numerically.  F.or  example,  if  we  have  avail- 

able  G^^(z)  * l^z)  * anc*  ^ » then  we  would  be 

able  to  obtain 


A* 


i,.(0  " *i.«  AVZ>  + Vl.m  A*l<*>  + ei+2.„  AV*>- 
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by  solvin8  «•«  3 equations  simultaneously.  In  the  case  at  hand, 
we  will  use  a site-diagonal  self-energy  since  all  we  have  at  our 
disposal  are  the  computed  G^(z;.  We  will  still  use 


x ik-(R  -R  ) 

G (2)  = — r e m 

N h - — _ 

K 2 - A*(z)  - e 
hence  l 


z-A*(z)  - e 


as  the  definition  of  the  off-diagonal  elements.  The  propagators 
decay  rapidly  with  distance  R^m  so  both  the  one,  two,  or  three  point 
curve  fitting  proceedures  will  probably  give  reasonably  equivalent 
results.  Now,  all  the  information  contained  in  our  previous  numeri- 
cal work  is  stored  in  A*(z)  , the  complex  proper  self-energy  part. 

It  is  of  interest  to  compare  Im  A*  with  Im  S in  order  to  see  how 

vi  A 

they  differ.  This  is  done  in  Figure  10  for  the  l-D  alloy  of  Figures 
5 and  6 . 

In  summary  we  have  presented  a relatively  simple  method  for  calculating 


the  eigenvalue  spectrum  of  a disordered  system,  one  that  avoids 
all  the  computational  pitfalls  of  self-consistent  methods.  This 
quasi-invariant  theory  is  not  only  highly  accurate,  but  also 
allows  the  bounds  on  the  frequency  spectrum  to  be  naturally 
determined  by  the  correlated  scattering  of  a local  group  . 

We  now  discuss  transport  and  develop  a formalism  that  allows 
our  numerical  output  to  be  used  in  approximations  that  con- 
serve particle  number  and  energy.  • 
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VI  Transport  tn  Disordered  Systems 


The  linear  response  of  the  current  to  the  electric  field  defines 
the  conductivity,  which  we  take  to  be  the  same  along  the  3 principal 

,14 

directions  in  our  simple  cubic  structure.  Following  Velicky  we  have 
in  our  single-particle  model 


a<0)  = ^ J*  dX 


where  e , the  electric  charge,  is  unity  and  f is  the  fermi  function. 
The  bracketed  term  is  short  hand  for 


<6(X-H>Pl6U-H)Pl>  - g <ol6(X-H)Pl6(X-H)Pl|a 


P.  is  the  momentum  operator  along  an  arbitrarily  chosen  principal  axis 
1 

and  the  long  bar  denotes  configuration  averaging.  Examination  of  Eq(i4) 
reveals  that  we  require  the  two-particle  correlation  function 


gLu(zi>22>*  <‘l  41  ferula  • 


We  can  relate  G to  G by  the  equation 


GLli(2l-22>-<:iJ(7l>C«i(:t2>+k5Gik<Zl)G<1t<22)-kr,MGp„Jr<ll-'2)-06) 


which  defines  the  index  structure  of  the  vertex  function.  This 


equation  is  exact  for  the  exact  G and  we  will,  use  it  to 

2~ 

define  G when  we  have  an  approximate  single  particle  green  function. 

We  can  place  restrictions  on  possible  vertex  functions  by 
requiring  the  conservation  of  charge  and  energy  in  the  presence  of 
a long-wavelength  disturbance.  This  leads  to  the  introduction  of  a 
new  operator 


K(VV  ' z^-H 


It  is  easy  to  show  that  K must  satisfy  the  following  Ward-type  identity: 

K(Z1,Z2)  = '2 l"-z  ' G(*2)]  (18) 

which  implies  a connection  between  the  one-body  operator  G and  the 

2 

two-body  green  function  G . Also,  one  can  show  with  the  above 
condition  that  the  linear  response  of  the  particle  number  and  energy 


to  a long  wavelength  disturbance  is  zero,  thus  ensuring  the  appropriate 
conservation  laws.  If  we  use  an  approximate  G,  then  we  must  construct  a 
K that  maintains  Eq . (18)  and  this  allows  us  to  relate  G to  the  vertex 
function  in  the  following  way:  We  let  *2  “*  z\>  then  Eq(18)  becomes 


Ktt(vV 


d ci<.(V 


The  configuration  averaged  resolvent  can  be  written  as 


■ <* 


z.  - T.  . -A  (z, ) 
l kin  op  1 


which  leads  to  the  equation 


(21) 


Vzi-‘l>  Vzl>  Vzi>  +kjj 


But 


dz. 


£ ^\q(zt}  r~dG.r(zl>1  £ dA*kq(zl> 


rp  dGpr<zl> 


dz. 


21  » Z1  ) 


rP  dGpr<zl>  ^ 1,1 


so  that 


Gij(zi)Gj/zi>+^  °ik<zi>Vzi>  d*Vzi>  Vzi’zi)-  ' 


kq 

rp 


dG  (z.) 
pr  1 


(22) 


In  G , we  let  m = j and  sum  over  all  j with  z^  = z^ : 


j GijjX(2l,Zl)  * Gij(zl)GjX(2l)  + ^ Gik(zl)Gqi(zl>  Skr,pq(zl»zl^Gpjjr(zl>zi: 
J J rp  J 

(23) 


Once  we  recognize  K.  (z. , z- ) = £ G... . (z.  ,z  ),we  find  that  the 

J.  1 j Xjj£  11’ 

vertex  must  satisfy  the  equation 


skr,pq<zl>zl>  - dn0<2l> 


dG  (zn) 
pr'  V 


(24) 


We  note  that  not  only  must  this  relationship  hold  for  the  exact 
vertex  function  and  self -energy  but  also  in  any  approximation  in 


in  which  it  is  desired  that  the  two-particle  correlation  function 

satisfy  the  Ward  Identity,  Eq.(18).  This  means  that  once  an  approximation 
* 

is  made  to  A , we  can  determine  transport  functions  that  allow  the 

conservation  laws  to  be  obeyed.  We  could  make  further  approxi- 
2 

mations  to  G but  we  would  then  have  no  guarantee  of  conserving 
charge,  energy,  etc.  The  Ward  Identity  is  useful  to  generate  a 
vertex  function  only  when  the  frequencies  are  the  same.  For 
the  case  at  handjA^Cz^)  = 6^q  A*(z^)  is  a function  of  the  site- 
diagonal  averaged  Green  function  so  the  vertex  is 

d A*(z1 ) 

Skr,pq(Zl»2l)  = 6kq6rp  dG"(Zl)  = 6kq6rp  H(Z1’Z1>‘  (25) 

For  f Z2»we  make  the  approximation 


Hkr,pq(zl’22)  " 6kq  6rp  (26) 

This  is  certainly  consistent  with  the  Ward  Identity, and  further- 
more, it  allows  us  to  show  that  contributions  to  the  conductivity 
from  the  vertex  correlations  then  vanish.  The  two-particle  correla- 
tion function  is  now 


GLji('l"z2)=Gij<zl>Gra/'2>+flk<.zl>Gkl(22>  S<V“2>*4.jr(V22>  • (27> 


To  find  the  conductivity  we  need, (6(\1-H)p16(X2“H)P2)  *or 


*12^1^2^  = ^ Pjm(m|6(\2“^)  U)  P«i(^l6(A^_H)  | j) 


ijjjm 


.... - — - - - 


(28) 


which  requires 


SjjemPjn>  P£i  GLj/XlA2)  " ijxJVPjfci  Gij(Xl)  G™£(X2) 


+.S  E p*i  Gik(XP  Gk£(X2^  -(xl»X2^  E Grrajr(Xl,X2^  Pjm*  (29) 
urn  k * r 


The  second  term  breaks  up  into 


E P£i  Gik(xP  Gk/X2^  E H(X1’X2^  Grmjr(Xl’X2)  p jm  = A * B * (30> 
i£k  * * jmr 


Let  us  transform  the  Wannier  sum  in  A to  a Bloch  sum.  Then  since 

/kip  |k')=mV  (k)6  , and  G i-s  diagonal  in  the  Bloch  repre- 

12  2 £k' 

sentation. 


A = m £ V0 (k)  G (X,)  G (x9) 

V*  2 kk  1 £k  1 


= m £ 


Be 


Sk2  (x1-A*(X1)-e_<)  (X2-A*0.2)-s 

k 


The  propagators  are  even  under  inversion  (k  -•  -k)  but  the  velocity  v2= 
is  odd, giving  us  zero, and  all  vertex  corrections  now  vanish. 

In  this  case,  the  fortunate  cancellation  of  vertex  corrections 
comes  about  as  a consequence  of  the  approximation  of  the  proper 
self-energy  by  a site-diagonal  quantity,  Eq.(12). 

Zero-Temperature  D.C.  Conductivity 

We  are  now  in  a position  to  evaluate  a(0).  Because  the 
vertex  corrections  vanish, 

I12(X,X)  = £ p)m  <ml6(\-H)U>  p*t  <il&(A-H)|j> 

ijXm 

= m2  E VL(k)  V2(l^)  <k(6(X-H)^>2  • 


(31) 

) 

k 


Be/Bk2 


(32) 


With  the  definition  <K 1 6 (X  - H) |k>  = -1  Ira  G (\)  , we  get  for 
_ _ " kk 

1=2  (o(ra)  = an(uO) 


o(0)  = 2 — dxr"^1]  E Vl(^)2  [ImG^Cx)]2  , 


where  we  have  included  a factor  of  two  for  the  two  possible  spin 

orientations.  At  T = 0,-  ilfXA)  = 6(X  “ p)  with  p = chemical  potential, 

dX 

and  the  conductivity  per  atom  is 


2 E Vn(k)2  [ IraG  (H)l2 

ti  N k 1 kit  ’ 


o(0)  -2_J^<se  [ JsiSkJ fijv.iE)2  i(E,  ).,3M 

TT  Cp-ReA*(p  )-E)2  + ImA*Cx+)2  Nk  1 V °5) 

This  natural  separation,  only  possible  for  a proper  self -energy 

independent  of  k >isolates  the  lifetime  and  energy  shifts  of  the 

single -particle  excitations  from  that  part  of  the  conductivity 

which  pertains  to  the  particular  lattice  under  study.  We  will 

concentrate  on  a 3-D  simple  cubic  lattice  with  e =~fcosk  -tcosk  4cosk  1 

r 3L  x y ZJ* 

-•221  2 * J 
Then,V^(k)  = sin  k^  = ~ (1  - cos  k^) . Consider  the  functions 

1 — g 
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v nik 

1 e x 

P»(E>  -iSrV 

k V 


for  which, we  have 


ij  5 VjtS)2  UE-eJ  = igj-  I»  [P2(E+>  - Pc(E+)3. 
k k \ 


4*  4* 

The  imaginary  part  of  both  P0(E  ) and  P (E  ) vanish  outside 

* o 

the  unperturbed  band  and 


r(0) 


= 1 Im 

9rr2 


J-l  dE  t 


IinA*(u+) 


(ii-ReA*(n+)-E)2+  ImA* (jj.  ' ) 


— ]2  [P2(E+)-Po(E+)].  (36) 


We  have  calculated  the  D.C.  conductivity  for  our  3-D 

alloy  in  order  to  illustrate  our  formal  results  . Generally, 
there  are  two  ways  in  which  the  D.C.  conductivity  can  vanish. 

If  the  density-of-states  at  the  Fermi  level  is  zero  then  so  is 
o(0)  . In  addition,  we  can  have  a finite  e (u) » but  a zero  mobility 
because  of  wave  function  localization.  Eq.(36)only  admits  a zero 
in  o(0)  if  p( is  zero  so  we  cannot  take  the  latter  possibility  into 
account  . The  conductivity  is  displayed  in  Figure  11  against  its 
respective  density-of-states.  There  is  a rather  direct  correlation 
between  the  magnitude  of  the  density-of-states  and  that  of  the 
conductivity.  This  relationship  is  understood  in 
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Figure  Captions  > 
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Li 


Fig.  I Real  and  Imaginary  parts  of  the  complex  self-energy  in  the  CPA. 

The  two  sets  of  curves  are  Indicative  of  the  type  of  results 
that  can  be  expected  from  this  approximation.  In  part  a,  we 

have  a situation  where  the  alloy  bands  overlap  while  in  b the 

* 

case  of  split  bands. 

Fig.  2 Comparison  of  1,3,  and  5 cluster  calculations  for  p(u>)  using 

|lm  I = *15  with  = + .5,  c = .5.  Background  j 

(histogram)  is  exact  results  of  Ref.  4. 

Fig.  3 Comparison  of  1,3,  and  5 cluster  calculations  for  p(u)  in  1-D 

• i 

using  |lm  I = *5  with  = + 1.0  and  c = .5.  Background 

(histogram)  is  exact  results  of  Ref.  4. 


Fig.  4 Curves  give  p(ui)  for  3 different  values  of  E tor  a 5 cluster 

tl  id! 

calculation  of  a l-D  alloy  with  = ± 1.0  and  c * 5.  The 
sharpest  peaks  are  associated  with  lowest  value  of  |lm  I? 

.4.  Other  values  are  .5  and  .6. 


Fig.  5 (a)  Density-of-states  for  a 50-50%  concentration  l-D  alloy  with 

V = 2.0,  V = -2.0.  Histogram  is  exact  calculations  from 

o D 

Ref.  4.  The  full  band  is  obtained  by  reflecting  the  portion 
shown  through  the  origin.  These  results  were  computed  from 
Im  £ * -.80  and  Re  £ obtained  from  Im  £ by  Eq .( 8), including  the 
exact  scattering  from  all  configurations  of  a central  atom  and 
its  4 nearest  neighbors.  (b)  Extension  of  the  above  results 
to  a cluster  of  7 atoms  using  the  same  £(u>).  This  result  is 
comparable  in  accuracy  and  wealth  of  detail  to  the  best  self- 
consistent  calculation  to  date,  Ref.  4. 


V 
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Figure  Captions 


Fig.  6 Curves  give  p(u>)  for  3 different  values  of  for  a 

5 cluster  calculation  of  a l-D  alloy  with  = + 2.0  and 

c = .5.  The  sharpest  peaks  are  associated  with  lowest 

value  of  llm  E , , = .4,  other  values  are  .8  and  1.0. 

1 trial  ' * 

Fig.  7 Comparison  of  7 cluster  calculations  of  p((u)  (dashed  line) 

using  a two  step  |lm  I (long  dashed  line),  with 

13 

numerical  work  of  Alben  £t.  al.  (solid  line)  who  solved 
the  Schrodinger  equation  for  an  8000  atom  3-D  tight  binding 
solid.  The  potentials  are  = + .75  with  impurity  concen- 
tration of  .1.  Small  horizontal  arrows  indicate  the  height 
to  which  their  peaks  rise. 

Fig.  8 Comparison  of  7 cluster  calculation  of  p ( ui)  (dashed  line) 
using  a 3 step  |lm  I (long  dashed  line)  with  results 


of  Alben 


et.  al . (solid  line).  V.  = + .75  and  c = .1 


for  this  3-D  tight  binding  alloy.  The  sharp  peak  in  Fig.  7 
at  U)  = .86  is  absent  because  of  the  coarser  energy  scale 
used.  Arrow  indicates  the  height  to  which  their  peak  rises. 

Fig.  9 Density-of-states  for  V = .75,  V = .75  in  a 3-D  simple 

A D 

cubic  lattice  with  concentration  c^  = .1.  Cowley  short 
range  order  parameter  or  is  -.07,  .3,  and  .7  respectively. 

|lm  E|  is  the  same  as  in  Fig.  7 (dot-dash  line)  and  ReS 
is  obtained  therefrom  by  use  of  Eq . (8). 
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Figure  Captions 


Fig.  10.  Solid  line  is  |lm  A (ou)  [ for  the  1-D  alloy  of  Figs.  5 and  6 
in  the  5 cluster  approximation  while  the  dashed  line  is  the 
corresponding  |lm  E («u)  |. 
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Fig-  11  Then  D.C.  conductivity  (dashed  line)  of  a 3-D  alloy  with 
Vj  = + .75  and  c = .1.  Density-of-states  taken  from 
Fig.  7 is  shown  in  the  solid  line. 
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